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Chapter  1 
Introduction 


1.1  General  Background 

The  recent  explosion  in  internet  traffic  proclaims  the  arrival  of  the  Information  Age  [1]. 
This  exploding  internet  traffic  volume  at  the  same  time  imposes  a  tremendous  challenge 
on  communication  engineers  to  meet  the  unlimited  demand  for  information  bandwidth, 
both  cable  based  and  wireless  [2].  On  the  other  hand,  this  demand  also  brings  us  oppor¬ 
tunities  that  we  were  not  able  to  imagine  before.  A  recent  prediction  asserts  that  “data” 
communication  will  soon  overtake  “voice”  communication  [3],  and  future  development 
is  heading  toward  personalization  and  mobilization  [4].  A  wireless  information  super¬ 
highway”  is  envisioned  to  fulfill  the  demand  for  multimedia  networks  with  multi-service 
requirements  [5-7]. 

At  present,  systems  of  wireless  communication,  mainly  cellular,  are  only  able  to  pro¬ 
vide  low-speed  data  communication  at  error  rates  that  are  far  from  acceptable  in  wireless 
network  connections.  This  is  because  these  systems  have  been  established  based  on  tech¬ 
nology  distant  from  theoretical  limits  [8].  One  of  the  main  reasons  for  this  drawback  in 
technology  is  the  hurdle  of  multipath  fading  channels  [9,10].  A  very  large  proportion  of 
the  data  bits  is  used  for  signal  integrity,  either  by  extensive  code  correction  or  multiple 
repetition.  The  cellular  system  we  use  today  is  designed  more  to  circumvent  the  problems 
imposed  by  the  multipath  fading  phenomenon  than  actually  to  solve  them.  Even  if  there 
are  efforts  to  overcome  the  problems  of  multipath-fading  channels,  those  techniques  applied 
are  relatively  primitive.  A  typical  example  is  to  repeatedly  transmit  a  message  through 
a  channel  so  that  a  correct  reception  of  the  message  can  be  determined  by  a  majority 
vote  (e.g.,  the  United  Kingdom  total  access  communications  system  repeats  its  message 
11  times  [11])-  Also,  to  avoid  excessively  rapid  fluctuation  of  the  received  signal  power 
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caused  by  fading,  a  specific  signaling  rate  is  chosen  according  to  the  rate  of  fading  so  that 
the  received  signal  can  be  regarded  as  approximately  constant  over  a  data  symbol  interval 
(i.e.,  the  slow  fading  assumption).  At  the  same  time,  the  signaling  rate  is  limited  by  severe 
problems  with  inter-symbol  interference  (ISI),  which  is  a  frequency  selective  problem  due 
to  multipath  transmission  [12].  However,  there  is  an  indication  that  high-speed  reliable 
communication  is  possible  if  the  user  is  willing  to  optimize  what  the  media  could  pro¬ 
vide  [13].  Because  of  this  belief,  considerable  research  effort  has  been  devoted  in  the  area 
of  multipath- fading  channel  communication  for  the  past  five  decades  [9,14-22]. 

Usually  multipath-fading  channels  are  merely  called  fading  channels  [23].  However, 
there  are  actually  two  basic  degradations,  i.e.,  multipath  and  deep  fades.  Multipath,  which 
results  from  reflections  of  the  transmitted  signal  from  reflective  surfaces  between  the  trans¬ 
mitter  and  the  receiver,  causes  delayed  and  scaled  versions  of  signals  to  be  superimposed 
at  the  receiver.  When  the  arrival  times  of  the  different  rays  are  of  the  same  order  of  mag¬ 
nitude  as  the  duration  of  the  transmitted  signals,  successive  signals  are  smeared  together, 
thus  resulting  in  ISI.  The  span  of  the  excess  delay  (i.e.,  the  time  between  the  first  and  last 
received  components  during  which  the  received  signal  power  falls  to  some  threshold  level 
below  that  of  the  strongest  component)  is  directly  related  to  the  physical  distortion  of  the 
signal  by  the  channel.  This  kind  of  distortion  varies  according  to  the  difference  in  frequency 
contain.  Frequency  selectiveness  is  the  term  we  use  when  the  distortion  rate,  the  channel 
gain,  and  the  phase  variation  change  significantly  for  a  small  variation  in  frequency.  When 
the  arrival  time  difference  is  comparable  to  the  period  of  the  carrier  frequency,  deep  fades 
result.  Deep  fades  are  the  phenomena  that  waves  of  different  phases  superimpose  and 
interfere  constructively  or  destructively.  This  kind  of  interference  may  cause  an  extremely 
low  received  signal  power.  Sometimes,  this  deviation  of  signal  power  can  be  more  than 
40  dB  [9]!  The  distance  between  nulls  is  approximately  0.33  meter  for  900  MHz  cellular 
systems.  Deep  fades  can  also  result  from  a  relative  motion  between  the  transmitter  and  the 
receiver,  which  also  results  in  superposition  interference  called  the  Doppler  effect  [24].  At 
very  high  frequency  (VHF)  and  ultra  high  frequency  (UHF),  a  vehicle  moving  at  50km/hr 
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or  14  m/sec,  will  pass  through  several  fades  in  a  second  [11].  Deep  fades  can  also  be 
caused  by  the  continuous  physical  changes  of  the  channel  media.  When  the  interference 
and  fluctuation  are  much  more  rapid  than  the  signaling  rate,  fast  fading  occurs. 

Frequently,  fading  and  multipath  phenomena  occur  simultaneously.  However,  there 
are  situations,  such  as  stationary  communications,  in  which  both  the  transmitter  and  the 
receiver  are  not  moving,  where  the  main  cause  of  degradation  of  communication  quality 
is  the  occurrence  of  multipath  transmission.  There  are  also  situations  like  communication 
between  high-speed  vehicles  in  an  environment  with  no  major  reflective  surfaces,  where  the 
main  cause  of  degradation  is  fading.  To  make  the  multipath-fading  channel  problem  more 
tractable,  it  is  helpful  to  consider  the  multipath  aspect  or  the  fading  aspect  separately, 
and  this  is  the  approach  we  use. 

1.2  Literature  Review 

The  history  of  fading  channel  studies  can  be  roughly  divided  into  three  periods.  The 
first  period  spans  from  the  early  1950’s  to  the  late  1960’s.  The  driving  force  behind 
this  research  was  mainly  the  impetus  for  the  development  of  long-distance  troposcatter 
communication  for  military  use  [14,15,25-27].  Tremendous  effort  was  invested  in  collecting 
real  data  to  illustrate  fading  channels  [28].  Some  practical  mathematical  models,  such 
as  uncorrelated  scattering,  were  introduced  [14, 17].  The  single  most  important  idea  for 
combating  fading  phenomenon,  “diversity,”  was  also  formulated  in  this  period  [29-36].  The 
beneficial  effect  of  diversity  vanishes  when  the  channel  has  a  very  high  correlation,  i.e.,  no 
fading  [22].  A  lot  of  work  was  also  performed  in  the  area  of  detection.  However,  systems 
were  generally  not  reliable  enough  for  some  computer  communication  applications.  This 
was  the  main  reason  why  coding  theory  was  introduced  [19]. 

The  second  period  spans  from  the  early  1970’s  to  the  mid  1980’s.  During  this  period, 
there  was  a  drop  in  interest  in  fading  channel  communication,  which  might  be  due  to 
the  high  demand  in  research  of  wire  and  optical  communications,  where  fading  is  not  an 
important  issue.  Even  so,  there  was  still  continuous  effort  in  the  information  and  coding 
theory  aspects  of  fading  channel  research  [37-40]. 
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The  third  period  spans  from  the  mid  1980’s  to  the  present.  During  this  period,  there 
has  been  a  huge  increase  in  scientific  activity  in  fading  channel  research.  Many  new  results 
have  been  reported.  The  main  driving  force  is  the  tremendous  increase  in  demand  for 
personal  mobile  wireless  communication,  of  which  the  following  references  specify  aspects 
of  the  research.  The  research  in  this  period  can  be  roughly  divided  into  three  categories, 
i.e.,  implementation,  characterization,  and  optimization. 

1.  Implementation: 

Many  new  implementation  schemes  have  been  derived  and  reported.  These  ap¬ 
proaches  can  be  further  subdivided  into  three  more  classes:  a)  Diversity,  b)  Detection 
and  Estimation,  and  c)  Combined  Methods. 

(a)  Diversity:  The  application  of  the  idea  of  diversity  exists  in  virtually  every  branch 
of  research  on  fading  channels,  either  implicitly  or  explicitly.  A  basic  motivation 
for  diversity  is  that  under  fairly  general  conditions,  a  channel  affected  by  fading 
can  be  transformed  into  an  additive- white-Gaussian-noise  (AWGN)  channel  by 
increasing  the  branches  of  diversity  [36,41-45].  There  are  mainly  three  types  of 
diversity,  i.e.,  time,  space,  and  polarization  diversity.  Time  diversity  is  usually 
achieved  by  using  some  kind  of  rearrangement  of  the  signaling  process  [46],  or 
by  means  of  coding  [35,41,47,48].  Space  diversity  can  be  achieved  by  using 
multiple  antennas  at  the  transmitter  and/or  the  receiver.  Polarization  diversity 
is  achieved  by  using  an  antenna  or  antennas  with  multiple  polarizations  [49-52]. 

(b)  Estimation  and  Detection:  There  are  hybrid  adoptions  of  classic  techniques 
of  estimation  and  detection  theory  for  fading  situations,  usually  under  specific 
settings  [53-56].  For  example,  there  are  applications  of  diversity  on  the  design 
of  detector  and  estimator  [36,57],  and  the  usage  of  robust  adaptive  algorithms 
in  equalizer  designs  for  fast  fading  channels  [58-66]. 

(c)  Combined  Methods:  These  methods  include  information  theory  and  spatial 
diversity  [67],  adaptive  code-division- multiple-access  (CDMA)  signaling  [68], 
space-time  modulation  [49],  and  the  applications  of  channel  coding  [69,70]. 
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2.  Characterization: 

Tremendous  effort  has  been  devoted  to  the  study  of  channel  characterization,  the 
polishing  of  necessary  mathematical  tools,  and  the  development  of  performance  eval¬ 
uation  techniques.  Furthermore,  theoretical  performance  bounds  have  been  derived, 
especially  for  matched  filter  detectors. 

(a)  Channel  Characterization:  Channel  characterization  is  generally  carried  out 
on  an  experimental  basis.  Measurement  of  transmitted  and  received  signals 
are  made  in  typical  environments  in  a  controlled  manner.  Statistical  models 
are  formulated  and  validated  against  the  measured  data  and  applied  to  system 
analysis  and  simulation.  Some  researchers  collect  their  own  data  and  construct 
their  own  channel  models  [71]  while  others  study  the  existing  experimental  data 
with  their  own  models  or  improve  previous  models  [72,73].  The  demand  for 
research  in  this  topic  has  also  inspired  the  development  in  specific  mathematical 
techniques  [46,74-76]  and  new  statistical  models  [77].  The  characteristics  and 
effects  of  the  channel  on  the  various  aspects  of  system  performance  (e.g.,  error 
rate  under  specific  settings,  channel  capacity,  etc.)  alone  have  also  induced 
much  interest  among  researchers  [78-81]. 

(b)  Simulation  Techniques:  In  the  design  of  a  communication  system,  a  system 
designer  may  want  to  ensure  that  the  performance  of  the  communication  link 
is  satisfactory.  To  ensure  the  system  is  performing  up  to  the  design  specifica¬ 
tions,  a  computer  or  hardware  simulation  is  required.  Therefore,  many  different 
simulation  schemes  have  been  proposed  for  different  settings  [76,82,83]. 

(c)  Performance  Analyses:  Performance  analyses  have  been  derived  for  uses  in 
nearly  all  common  fading  situations  and  standard  communication  techniques. 
These  analyses  are  often  case-specific  [42,43,45,46,51,74,84-92]. 

(d)  Performance  Bounds:  Different  performance  bounds  that  have  been  derived 
now  set  guidelines  to  obtain  optimal  performance  for  different  fading  channels 
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[93,94].  Among  different  performance  bounds,  the  matched-filter-bound  (MFB) 
has  received  special  interest  from  researchers  [34,95-99].  The  MFB  characterizes 
the  detection  capability  of  a  single  communication  pulse  transmitted  in  isolation 
so  that  potential  effects  of  ISI  can  be  negligible.  Although  the  MFB  is  a  lower 
bound  over  a  communication  channel  that  may  not  be  practically  realizable, 
it  is  a  common  performance  measurement  for  communication  systems  [96,97]. 
There  are  reasons  for  saying  MFB  is  an  optimal  performance  bound.  First, 
it  is  assumed  that  all  the  information  about  the  channel  is  known,  which  is 
not  feasible  for  a  real  fading  channel.  Second,  transmitted  pulses  are  separated 
sufficiently  so  that  no  ISI  occurs.  In  other  words,  the  error  rate  is  determined 
by  assuming  only  one  data  pulse  is  transmitted,  which  is  also  not  physically 
realistic. 

3.  Optimization: 

A  great  amount  of  valuable  work  has  gone  into  assessing  the  information-theoretic 
limits  of  fading  channels.  Tremendous  effort  has  been  devoted  to  the  derivation  of 
channel  capacity  and  the  search  for  optimal  signals  to  achieve  its  optimal  capacity. 

(a)  Channel  Capacity:  The  capacity  with  or  without  the  knowledge  of  channel  state 
information  has  been  studied  [67,100,101],  or  with  some  further  constraints,  such 
as  the  delay  limited  to  real  time  voice  and  video  communication  [102].  Gold¬ 
smith  [100]  and  Caire  [101]  interpreted  the  optimal  power  adaptation  as  “water¬ 
pouring”  in  time  or  frequency,  while  for  multipath-fading  channels,  Telatar  [103] 
interpreted  these  as  “peaky”  in  time  or  frequency.  For  multi-access  usage  of  the 
fading  channel,  there  is  a  solution  analogous  to  the  “water-pouring”  interpreta¬ 
tion  for  the  single-user  case  to  maximize  the  overall  capacity  [104].  In  general, 
the  capacity  as  well  as  the  capacity-achieving  distribution  implies  that  there  are 
some  underlying  structures  of  optimal  coding/signaling  [105]. 

(b)  Optimal  Power  Allocation:  With  the  assumption  of  block-fading  additive- white- 
Gaussian-noise  (BF-AWGN)  and  the  knowledge  of  channel  states,  an  optimal 
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power  control  for  both  the  transmitter  and  the  receiver  is  derived  with  or  with¬ 
out  the  presence  of  a  transmission-delay  constraint  [106, 107].  A  discrete  ob¬ 
servable,  usually  orthogonal  projections  of  the  tested  signal  onto  a  family  of 
specific  functions,  has  been  used  to  reduce  the  receiver  complexity  in  the  study 
of  the  performance  for  difference  classes  of  signaling  [108].  Also,  there  are  more 
general  results  for  signaling  properties  for  fading  channels  [109, 110]. 

1.3  Methodology 

The  research  in  this  dissertation  utilizes  the  idea  of  separating  the  multipath  aspect 
from  the  fading  aspect  of  a  fading  channel  in  order  to  reduce  the  model  complexity.  The 
goal  of  this  dissertation  is  to  characterize  the  optimal  signal  and  detector  of  a  fading  chan¬ 
nel,  which  may  be  either  fast  fading  or  slow  fading,  when  its  channel  statistical  properties 

are  known. 

The  word  “optimal”  used  in  this  dissertation  denotes  maximized  signal-to-noise-ratio 
(SNR).  Improved  SNR  is  related  to  but  not  necessarily  equivalent  to  improved  error 
rate  [92].  Optimization  of  SNR  does  not  necessarily  mean  optimization  of  the  error  rate. 
Optimizing  in  terms  of  error  rate  is  the  ultimate  goal  for  all  communicators.  It  involves 
sophisticated  modeling  of  the  entire  communication  system  including  the  optimal  detec¬ 
tor  corresponding  to  the  autocorrelation  function  (which  is  a  stochastic  process  itself)  of 
the  fading  channel.  At  the  present  stage  of  research,  we  do  not  try  to  tackle  this  prob¬ 
lem,  taking  on  instead  the  more  tractable  optimization  of  SNR.  This  is  analogous  to  the 
SNR-maximizing  derivation  of  the  conventional  matched  filter,  which  also  happens  to  re¬ 
sult  in  minimal  bit  error  rate  in  binary-phase-shift-keying  (BPSK)  systems.  In  any  event, 
improvement  in  received  SNR  will  result  in  an  improvement  in  the  bit  error  rate. 

We  first  characterize  optimal  detectors,  then  study  optimal  signaling.  Because  of  our 
model  setting,  these  two  optimization  problems  are  essentially  the  same  mathematically, 
which  significantly  reduces  the  difficulty  of  the  derivation  and  understanding  of  the  char¬ 
acterization  of  the  optimal-signal-filter-pair  for  fading  channels.  For  the  sake  of  clarity  of 
later  discussion,  we  call  this  optimal-signal- filter-pair  a  consonant  pair  (OP). 
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We  use  functional  analysis  techniques  to  derive  the  characteristics  of  the  optimal 
detector  and  signal  in  £2-1  We  then  further  analyze  the  properties  of  the  optimal  filters 
and  signals,  and  present  some  algorithms  that  help  us  to  find  these  properties. 

Without  the  assumption  of  slow  fading,  the  matched  filter  is  no  longer  necessarily 
the  optimal  detector  [111,  112].  As  a  result,  the  SNR  of  the  received  signal  will  depend 
strongly  on  the  channel  statistical  properties,  the  corresponding  transmitted  function,  and 
the  detection  filter.  Therefore,  with  proper  functional  analysis  techniques,  we  can  derive 
some  basic  properties  of  the  signal  and  filter  pairs  to  maximize  the  SNR. 

There  are  two  main  results  from  this  research.  First,  maximizing  the  SNR  of  the 
detected  signal  is  a  key  factor  to  lowering  the  error  rate  in  communication,  which  implies 
a  higher  possible  data  rate,  increased  channel  re-usability,  smaller  appliance,  etc.  Second, 
a  better  characterization  of  the  optimal  signal  enables  us  to  have  better  resource  manage¬ 
ment,  e.g.,  if  we  know  that  the  support  of  the  optimal  signal  is  related  to  the  support  of 
the  autocorrelation  function  of  the  channel,  we  will  not  slow  down  our  transmission  rate 

by  using  long  signals. 

1.4  Contributions  to  the  Field  of  Knowledge 

After  an  intensive  literature  search,  we  believe  the  following  contributions  to  the  field 
of  optimal  detector  and  signaling  in  fading  channels  are  original  to  this  dissertation: 

1.  A  characterization  of  the  optimal  detector  for  known  signals  in  fading  channels.  We 
find  that  the  well  known  matched  filter  and  RAKE  receiver  are  limiting  cases  of  our 

result. 

2.  A  characterization  of  the  CP  in  terms  of  maximum  SNR  for  given  fading  channels. 
Not  only  is  fading  (multiplicative  and  multipath)  considered  as  a  channel  effect,  but 
also  the  effect  of  synchronization  uncertainty  on  the  optimal  signal/detector  design. 


i£a(R)  is  the  collection  of  all  Lebesgue  measurable  functions  f  :  R  ->  C,  for  which  |/|2  is  Lebesgue 
integrate,  i.e.,  f°°  \f\2  <  oo.  In  this  dissertation,  we  use  C2  to  represent  £2(R). 

J  —  OO 
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3.  An  iterative  algorithm  for  finding  CPs  that  can  be  used  for  communication  through 
fading  channels.  With  land-mobile-fading  channels,  the  signals  obtained  with  our 
procedure  are  demonstrated  to  have  a  performance  significantly  better  than  that 
of  the  traditional  flat-top  pulse  and  raised  cosine  pulse  in  a  PPM  pair  in  terms  of 
probability  of  error  rate  at  a  comparable  signaling  rate. 

4.  A  closed-form  expression  for  the  probability  density  function  and  the  cumulative 
distribution  function  of  quadratic  form,  g  =  vHQv,  where  v  is  a  complex  Gaussian 
vector  with  E(v)  =  0,  and  Q  is  a  Hermitian  matrix.  The  probability  density  function 
and  the  cumulative  distribution  function  are  expressed  in  terms  of  the  eigenvalues  of 
MQ,  where  M  is  the  covariance  matrix  of  v.  The  closed-form  expression  is  applicable 
to  many  general  cases.  In  this  dissertation,  we  apply  the  equation  to  the  calculation  of 
the  probability  of  error  rate  for  quadratic  receivers.  As  a  corollary  to  this  analysis,  an 
interesting  partition  of  unity  by  sequences  is  presented.  Prom  the  above  cumulative 
distribution  function  of  the  quadratic  form  of  normal  vectors,  we  are  also  able  to 
derive  an  expression  of  partition  of  unity  which  is  shown  to  be  very  useful. 

The  results  of  this  dissertation  are  primarily  analytical,  with  the  structure  of  the  CPs 
being  developed.  Nevertheless,  the  concepts  derived  herein  have  been  implemented,  and 
a  description  of  this  implementation  is  provided.  Also  a  modest  set  of  computed  results 
is  presented  to  demonstrate  the  nature  of  the  solutions  and  the  algorithms.  While  an 
exhibitive  demonstration  is  not  provided,  the  tools  described  here  can  be  used  in  a  variety 
of  circumstances. 

1.5  Dissertation  Structure 

In  chapter  2,  we  illustrate  the  fading  channel  model  that  is  commonly  used.  Following 
this,  we  introduce  the  model  which  we  have  used  in  our  research.  In  chapter  3,  we  discuss 
the  optimal  detection  of  a  given  signal  in  fading  channels.  Then,  we  relate  the  optimal 
signal  to  its  optimal  detection.  In  chapter  4,  we  develop  some  algorithms  that  we  can  use  to 
generate  CPs.  As  we  only  characterize  a  single  CP,  while  a  realistic  communication  requires 
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at  least  two  different  signals  to  represent  “0”  and  “1” ,  in  the  last  section  of  chapter  4,  we 
discuss  a  procedure  to  find  two  conjugate-consonant  pairs  which  can  be  used  for  binary 
transmission.  In  chapter  5,  we  derive  a  method  to  calculate  the  probability  of  error  for  our 
channel  model  and  compare  the  CPs  obtained  in  our  research  to  other  traditional  signals. 
Next,  we  discuss  the  issues  of  implementing  our  algorithms  with  discrete  samples  in  chapter 
6.  Results  of  the  implementation  are  presented  in  chapter  7  along  with  observations  and 
discussions.  We  conclude  this  dissertation  in  chapter  8  with  a  list  of  topics  for  further 


research. 
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Chapter  2 
Channel  Models 

2.1  Introduction 

Fading  channels  were  first  modeled  in  the  1950s  and  1960s  as  a  large  number  of 
“scatterers”  located  at  random  points  within  the  propagation  path  [16, 18].  This  idea  was 
primarily  applied  to  over-the-horizon  (troposcatter)  communications  covering  a  wide  range 
of  frequency  bands.  To  reach  beyond  the  horizon,  clouds  of  particles  in  the  troposphere 
were  used  as  reflectors  for  the  radio  waves.  The  differences  in  path  length  between  the 
large  number  of  scattered  waves  give  rise  to  Rayleigh  fading  if  there  is  no  dominant  direct 
component  [27].  Mobile  wireless  systems  of  recent  interest,  which  are  mainly  local  and  with 
finite  number  of  reflection  paths,  experience  fading  effects  that  are  somewhat  different 
than  those  mentioned  above.  However,  these  early  models  are  still  quite  useful  to  help 
characterize  fading  effects  that  we  are  facing  nowadays. 

The  groundwork  of  modeling  fading  channels  was  mathematically  laid  out  by  Price  in 
the  mid  1950’s  [14,15].  For  the  following  half  century,  there  were  tremendous  contributions 
to  this  field.  In  a  recent  survey  paper  written  by  Biglieri,  Proakis,  and  Shamai,  549 
citations  reporting  the  state-of-the-art  achievements  in  research  on  fading  channels  were 
included  [22].  Among  those  that  are  not  included  in  their  citations  is  a  paper  written 
by  Bello  who  introduced  a  simple  way  to  model  the  fading  phenomenon  with  the  notion 
of  wide-sense  stationary  uncorrelated  scattering  (WSSUS)  [17].  Uncorrelated  scattering 
means  the  attenuation  and  phase  shift  of  the  channel  associated  with  different  path  delays 
are  not  correlated. 

2.2  Fading  Channel  Model 

A  widely  accepted  mathematical  modeling  of  the  fading  channel  can  be  found  in 
Proakis’s  book  [21].  A  fading  channel  is  viewed  as  a  continuous  time-varying  filter  with  a 
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baseband-equivalent  impulse  response  of  c(r;t),  where  c(.; .)  is  a  complex  valued  function, 
r  is  the  delay  in  response  and  t  is  the  time  variable.  If  a  signal  s(t)  is  sent  through  this 
channel  with  AWGN,  we  will  have  at  the  receiver  side,1 

r(t)  =  J  c(T\t)s(t -r)d,T +  n(t),  (2.1) 

where  r(i)  is  the  received  signal  and  n(t)  is  the  AGWN. 

Assuming  c(r;  t)  is  wide-sense  stationary,  the  autocorrelation  function  of  c(r;  t)  can 

be  defined  as 

MW,  At)  =  ^[c*(ri;t)c(T2;i  +  At)].  (2.2) 

With  uncorrelated  scattering,  we  have 

4>c(ti;  At)S{n  -r2)  =  ^E[c*(ri;t)c(r2;t  +  At)].  (2.3) 


Setting  At  =  0,  the  expression  0c  ( Ar;0)  in  (2.3)  is  usually  called  the  Multipath  Inten¬ 
sity  profile  or  the  Delay  Power  Spectrum  of  the  fading  channel.  We  further  define  Tm 
to  be  the  length  of  time  that  0c  (At;  0)  is  essentially  nonzero,  and  call  it  the  Multipath 
Spread  of  the  fading  channel. 

Taking  the  Fourier  transform  of  the  channel  impulse  response,  we  have 

C(f‘,  t)  =  J  c(r;t)e~^TdT.  (2.4) 


The  corresponding  autocorrelation  function,  with  WSSUS,  is 
<M/i,/2;  At)  =  \E[C*(h-,t)C(h-,t  +  At)) 

=  J  ^[c*(ri;t)c(r2;t  +  At)]ej2^lTl_^2T2Widr2 

=  //  4>c{n’,At)S{n  -T2)^{hn~hTi)dTldT2 

=  I  Mn;A  t)e^A-AfrdrL 
=  j  <Mti;a i)e-;2*A/ndTi 


(2.5) 


=  0c  (A/;  At) 


throughout  this  work,  integrals  stated  without  limits  are  assumed  to  be  over  the  interval  (-00,00). 


Setting  A t  =  0,  we  obtain  (j>c(Af;  0);  this  is  the  autocorrelation  function  in  the  frequency 
variable.  We  define  the  length  of  the  support  of  (j>c(Af;0)  that  is  essentially  nonzero  as 
(A f)c-  Usually,  (A f)c  is  denoted  as  the  Coherence  Bandwidth,  which  is  a  measure  of 
the  frequency  coherence  of  the  fading  channel.  To  illustrate  this  idea,  we  would  imagine 
a  signal  with  a  bandwidth  less  than  (A f)c  being  transmitted  through  a  fading  channel. 
Because  all  frequency  components  will  be  affected  essentially  the  same  by  the  channel, 
this  channel  is  said  to  be  frequency-nonselective.  However,  if  another  signal  with  a  band¬ 
width  larger  than  (A f)c  is  being  transmitted  through  the  same  fading  channel,  different 
frequency  components  will  be  affected  in  different  ways  by  the  channel.  Then,  the  channel 
is  said  to  be  frequency-selective.  In  this  case,  the  signal  can  be  severely  distorted  by  the 
channel.  Another  useful  result  due  to  the  fact  that  ^c(Ar;0)  and  <f>c(Af;0)  are  Fourier 
transform  pairs  is  that 

(A f)c  «  (2-6) 

J-m 

To  study  the  time  variation  of  the  channel,  we  return  to  0c(A/;  At)  defined  in  (2.5). 
By  setting  A /  =  0,  we  obtain  the  autocorrelation  function  of  the  channel  for  each  different 
frequency  component  by  time  averaging.  Now,  the  length  of  the  support  of  (f>c{ 0;  At)  that 
is  essentially  nonzero  is  defined  as  (A t)c,  which  is  called  the  Coherence  Time  of  the 
fading  channel.  So,  in  a  digital  communication  system,  if  the  signaling  period  is  shorter 
than  (At)c,  we  may  assume  that  the  channel  is  essentially  constant  for  individual  signals, 
and  the  situation  of  slow  fading  occurs.  On  the  other  hand,  if  the  signaling  period  is  longer 
than  (At)c,  we  may  no  longer  assume  that  the  channel  is  constant  for  individual  signals, 
and  fast  fading  occurs. 

In  practice,  (Af)c  is  very  difficult  to  measure  directly.  So,  a  related  function  is 
employed.  This  new  function  is  called  the  scattering  function  of  the  fading  channel.  It  is 
defined  as 

S(t;  0  =  J  1 4>c(Af-,At)e-2^e~2^AtdAfdAt,  (2.7) 

where  r  is  the  time  delay,  and  (  is  the  Doppler  frequency  as  the  channel  varies  with  time 
(caused  by  physical  motions  of  the  media,  the  transmitter,  and  the  receiver,  etc).  The 
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scattering  function  of  the  channel  provides  us  with  a  measure  of  the  average  power  output 
of  the  channel.  When  r  =  0,  the  length  of  the  support  of  5(0;  ()  that  is  essentially  nonzero 
is  called  the  Doppler  spread  Bd  of  the  fading  channel.  We  have  the  following  relationship 
between  the  channel  coherence  time  and  the  Doppler  spread 

(A t)c  «  Yd'  (2’8) 

The  time  variation  autocorrelation  function,  </>c( 0;  At),  can  be  obtained  from  Sc{()  = 
5(0;  C),  and  is  called  the  Doppler  power  spectrum  of  the  fading  channel.  Some  of  the 
common  spectra  and  autocorrelation  functions  are  listed  in  Table  2.1  [9,84,92, 112—114]. 


Table  2.1:  Spectra  and  Autocorrelation  Function  of  Fading  Process. 


Denotation 

Spectrum  Sc  (0 

Autocorrelation  Function  Rc{t ) 

1.  Rectangular 

l£l  <  Bd, 

sin  2-k  Bdr 

27T  BdT 

2.  Gaussian 

Kexp{jjp)/y/ifBd 

K  exp[-(nBdT)2] 

3.  Land  Mobile 

K 

7T(C2-Bd2)l/2 

KJ0(2irBdT) 

4.  First-Order  Butterworth 

K 

nBd(l  +  t2/Bd2) 

K  exp[— 27rBd|r|] 

5.  Second-Order  Butterworth 

K 

1  +  16b7 

K  exp[  —^} 

x  (cos  ^1  +  sin  JfijP) 

Note:  Jo(-)  is  the  zero-order  Bessel  function  of  the  first  kind,  and  if  is  a  constant. 


2.3  Categories  of  Degradation 

Different  fading  phenomena  introduce  different  degradations  in  communication  sys¬ 
tems.  These  degradations  are  summarized  in  Table  2.2  [23]. 
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Table  2.2:  Categories  of  1 

Degradation  Accordinj 

to  Different  Fading  Phenomena. 

Frequency-nonselective 

Frequency-selective 

Slow  fading 

Fast  fading 

loss  in  SNR 

ISI  distortion, 
pulse  dispersion, 
irreducible  BER 

low  Doppler, 
loss  in  SNR 

high  Doppler, 
PLL  failure, 
irreducible  BER 

2.4  Model  for  Erequency-Nonselective  Fading  Channels 

Let  us  revisit  (2.1),  i.e., 


r{t)  = 


j c(T-,t)s(t  -  T)dr  +  n(t). 


By  Parseval’s  relation,  we  have 

r(t)  =  J  C(/;i)5(/)^td/  +  n(t).  (2.9) 

Assuming  frequency  nonselectiveness,  we  have  C(/;i)  =  C(0;  t),  and 

r(t)  =  J  C(/;i)5(/)e^d/  +  n(t) 

=  C{ 0;  t)  j  S(f)e>2vftdf  +  n(t)  (2‘10) 

=  C{0;t)s(t)  +  n(t). 

The  transfer  function  C(0;t)  for  a  frequency-nonselective  fading  channel  may  be  ex¬ 
pressed  as 

C(0;t)=a(t)=o(t)e"^W,  (2.11) 

where  the  random  process  a(t)  represents  the  envelope,  and  the  random  process  6{t)  rep¬ 
resents  the  phase  of  the  the  transfer  function.  In  this  case,  the  channel  is  also  called  a 
multiplicative  fading  channel.  Note  that  a  multiplicative  channel  can  be  both  fast-  and 
slow-fading.  For  different  fading  channels,  a(t)  can  acquire  different  statistical  distribu¬ 
tions,  from  which  each  different  fading  channels  obtains  its  name.  Among  common  distri¬ 
butions  for  a(t)  are  the  Rayleigh  distribution,  the  Rician  distribution,  and  the  Nakagami 
distributions  [92]. 
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2.5  Model  Used  in  This  Dissertation 

The  approach  that  is  employed  in  this  dissertation  for  modeling  the  general  fading 
channel  is  first  to  assume  frequency-nonselectiveness,  then  to  introduce  countably  many 
multiple  paths  when  we  want  to  study  the  frequency-selectiveness  aspect  of  the  fading 
channel.  This  approach  agrees  with  the  assumption  of  independent  fading  paths  put  for¬ 
ward  by  Bello  [17].  Furthermore,  this  approach  provides  a  model  that  is  general  enough 
to  represent  most  “real”  fading  channels  that  are  used  in  practice.  At  the  same  time,  the 
model  also  allows  us  to  separate  the  fading  and  multipath  aspects  of  a  fading  channel  and 
deal  with  them  separately. 

Fig.  2.1  depicts  the  system  that  we  are  studying.  It  is  a  generic  communication  system 
for  fading  channels  with  a  correlation  detector. 

When  we  assume  frequency-nonselectiveness,  a  fading  channel  becomes  simply  a  mul¬ 
tiplicative  Gaussian  channel,  i.e.,  when  s(t)  is  sent,  we  receive 

r(t )  =  a(t)s(t )  +  n(t).  (2.12) 

Without  loss  of  generality,  we  may  assume  a(t)  has  a  unit  variance,  i.e.,  E[a(t)a*(t)]  - 
E[a(t)]E[a*  (t)]  =  1.  For  a  general  multipath  fading  channel,  we  receive 

r(t)  =  ^2 &“*(*)«(*  “  &)  +  (2-13) 


_ 

transmitter 


channel 


receiver 


Fig.  2.1:  System  model  of  a  fading  channel. 
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where  ai(t)  has  an  unit  variance,  &  is  the  path  gain  for  path  i,  and  fa  is  its  path  delay. 
The  autocorrelation  functions  of  a.i{t)  =  Pi(t)  +  where  both  Pi(t)  and  qi(t)  are  real 

random  processes,  is 

Rai{r,  A)  =E[ai(T)a*(  A)] 

=E[(pi(T)pi( A)  +  qi(T)qi( A))  -  j(pi(r)<?i( A)  -  gi(r)p;(A ))]  (2.14) 

=E[{pi{T)pi(  A)]  +  E[qi(T)qi(  A)] 

where  it  is  assumed  that  the  real  and  imaginary  parts  are  independent.  The  autocorrelation 
Rai  (r,  A)  is  real  and  symmetric  about  r  and  A. 

When  ai(t)  is  Rayleigh,  i.e.,  pi{t)  and  qi{t)  are  zero  mean  Gaussian  processes,  and 
when  we  further  consider  that  the  different  fading  paths  are  only  caused  by  local  reflections, 
it  is  safe  to  say  that  the  ai(t)  all  have  the  same  normalized  statistical  characteristics,  i.e., 
we  have  for  some  autocorrelation  function  f?Q(r,  A)  that  Rai (r,  A)  —  Ra(r,  A)  for  all  i.  The 
above  simplification  cannot  be  applied  in  general  to  other  fading  channels.  For  example, 
if  ai(t )  is  Rician  distributed,  the  directions  of  the  dominant  components  are  different  for 
each  of  the  multiple  paths,  and  hence  the  autocorrelation  functions  may  be  different.  To 
simplify  the  derivation,  we  will  assume  a(t)  to  be  Rayleigh  in  this  dissertation.  However, 
our  discussion  is  readily  extendible  to  other  fading  channels. 

If  we  do  not  assume  perfect  synchronization  when  a  known  signal  so(t)  is  transmitted, 
the  received  signal  at  the  output  of  the  detector  filter  is 

r  =  f  £iai(t)s0(t  -  fa-  p)  +  n(t)]h{t)  dt, 

J  i 

where  p  is  a  random  delay  with  density  function  p{p).  We  assume  throughout  this  work 
that  h(t)  has  unit  energy,  i.e.,  f  h(t)h*(t)  dt  =  1. 


18 


Chapter  3 

Optimal  Detection  and  Signaling 

Up  to  the  present,  research  about  fading  channels  has  been  performed  mostly  under 
the  slow  fading  assumption,  i.e.,  the  channel  is  considered  to  be  constant  for  up  to  a  few 
symbol  times.  Under  this  assumption,  a  fading  channel  behaves  as  an  AWGN  channel 
for  each  symbol.  It  is  also  generally  accepted  that  matched  filter  detectors  produce  good 
performance.  However,  when  the  slow  fading  assumption  is  lifted,  the  idea  of  optimal  de¬ 
tection  becomes  more  complicated.  Fast  fading  causes  variations  within  a  symbol  interval, 
so  that  the  model  applied  in  slow  fading  becomes  inaccurate,  and  conventional  matched 
filter  detectors  become  suboptimal.  In  this  chapter,  we  characterize  optimal  detection  and 
signaling  without  the  assumption  of  slow  fading. 

3.1  Optimal  Detection 

For  a  single-path  frequency-nonselective  fading  channel  with  random  amplitude  func¬ 
tion,  a{t),  when  a  known  signal  so{t)  is  transmitted  in  a  Gaussian  channel,  the  received 
signal  is 

r(t)  =  a{t)so{t  -  p)  +n(t),  (3.1) 

where  p  is  a  random  delay  with  probability  density  function  p(p)  due  to  uncertainty  in 
the  synchronization,  and  n(t )  is  AWGN  with  a  power-spectral-density  (PSD)  No/ 2.  For  a 
multipath  fading  channel,  the  received  signal  is 

r{t)  =  &a»(*)a  o(*  -&-/*)  +  n(t),  (3.2) 

i 

where  &  is  the  channel  gain,  and  fc  is  the  path  delay  for  path  i.  We  would  like  to  find  the 
optimal  h(t),  with  unit  energy  /  h(t)h*(t)  dt  —  1,  that  maximizes  the  received  SNR. 
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3.2  Characterization  of  Optimal  Detectors 

The  received  signal  at  the  output  of  a  linear  detection  filter  can  be  presented  as 

r  =  J[Y.  £iOti(t)so(t  -  fa  -  p)  +  n(t)]h(t)  dt ,  (3.3) 

with  the  corresponding  signal  and  noise  components  at  the  output  of  the  detection  filter 


being 


=  y  dt>  and  rn  =  j  n{t)h{t)  dt.  (3.4) 


Since  h(t)  has  unit  energy,  we  have 


The  variance  for  the  signal  part  is 


2n  _  No 


E(\rn \2)  = 


E{\rs\2)  =E[J J 

=  m)so(A  fj)]h{T)h*{\)  dr  dX. 


Since  a(t)  and  s0(t)  are  independent,  we  have 


E(\rs\2)  =  / J £i£jE[oii{T)a*-(\)]E[so(T  —  fa  —  p)sq(X  —  fa  ~  m)]Mt)^*(A) 

i  j 

(3.7) 

Because  of  the  assumption  that  different  paths  are  independent,  we  obtain 
£(|rs|2)  =  Y  J  J  \Zi\2E[ai(T)a*{\)}E[so(T-fa-n)so{\-fa-n)]h(T)h*{\)  dr  dX.  (3.8) 


Defining 


®«(t,  A)  =  E[sq{t  -fa-  p)s*0(X  -fa-  p )] 


=  J  so(t- fa-  ix)s*Q{  X  -fa-  p)p(p)  dp, 


(3.8)  can  be  rewritten  as 


£(N2)  =  ?//  |&|2f?a(r,  A)tf i(T,  A )h(r)h* (A)  dr  dX 

=  JJ  Ra(r,X){Y\^\2Mr,m(T)h*(X)drdX. 


(3.10) 
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Let 


^(r,A)  =  x;i6|2^(T,A). 

i 

Then,  the  variance  for  the  signal  part  may  be  presented  as  a  functional  of  h. 

E(\rs\2)  =  //  Ra{r,  A)\?(t,  X)h(r)h* (X)  dr  dX 


(3.11) 


(3.12) 


=  F(h). 


Note  that  both  T(r,  A)  and  Ra(r,  A)  are  symmetric  in  r  and  A. 

By  applying  the  method  of  Lagrange  multipliers  to  maximize  F{h )  under  the  con¬ 
straint  of  H{h)  =  f  h(t)h*(t)  dt  =  1,  we  have  the  following  Lagrangian  function 

J(h)  =  F(h)-'yH{h),  (3.13) 

where  7  is  the  Lagrange  multiplier.  To  find  the  extrema,  we  take  the  derivative  with 
respect  to  h  (see  Appendix  A)  of  the  Lagrangian  and  set  it  to  zero,  i.e., 


>Afw  =  7jL*w, 


(3.14) 


where  h0(t)  is  an  extremum  of  F{h).  Prom  the  example  in  Appendix  A,  we  know  that 

dSM( r).« 

ah 

Uo{t)  +jv0{t), 


vv  uv^i.  - - —  -  —  /  a  *  * 

dg(/i0)(r)  _  NoW)  we  try  t0  ^°)(r).  For  eSi,  and  V  €  £2,  and  h0(<)  = 

dh  dh 


(3.15) 


.F((ito  +  eV)  +jv 0)  -  jP(/io) 

=  JJ  Ra(T,X)-${T,X){{uo  +  eV)+jvo)(T){{u0  +  eV)-jvo){X)dTdX-F{ho) 

=  JJ -Rq(t,  A^r,  A){(uo  +  eV)(r)(ito  +  el^)(A)  -  j{uo  +  eF)(r)uo(^) 

+  jv0(T)(uo  +  eV)(A)  +'y0(r)vo(^)}  dr  dX  -  F(h0). 

By  neglecting  those  terms  containing  e2,  (3.15)  becomes 
F((uq  -t-  eV )  +  jvo)  -  F(ho) 

«  JJ  i?Q(r,A)^(r,  A){uo(r)uo(A)  +  eF(r)uo(A)  +  euo(T)F(A) -^(^^(A) 

-jeV{r)v o(A)  +  jvo(r)u0(A)  +  jev0(T)V(X)  +  v0(t)vo(X)}  dr  dX  -  F(/io) 

=  JJ  j Rq(t,  A)\H(t,  A){eF(r)u0(A)  +  €U0(t)F(A)  -jeV(r)v0(A)  +  jeu0(r)V(A)}  dr  dX. 


(3.16) 
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As  both  Ra  and  are  symmetric  in  r  and  A,  (3.16)  may  be  further  simplified  as 
F((«o  +  eV)  +  jv o)  -  F(ho) 

=ej I  Ra(T,Xmr,X){uo(X)+uo(X)-jvo(X)+jvo(X)}V(r)dTdX  (3.17) 

=2e  J J  Ra(T,X)-${T,X)u0{X)V{T)  dr  dX. 

From  Appendix  A,  we  have 

aF(M(r)  _  2  f  Ra(r,  X )¥(r,  A)u0(A)  dA.  (3.18) 

du  J 


Similarly, 


(3.19) 


(3.20) 


(3.21). 


^iM(r)  =2  f  i?a(T,A)tf(r,A)vo(A)  dA. 
dv  J 

From  (3.18)  and  (3.19),  we  have 

=  f  RJ,T,X>1(T,\)h-Q(\)  ik. 

So,  from  (3.14), 

7 /ij(r)  =  j  Ra{r,X)^{T,X)h*0{X)  dX. 

Also, 

F(|rs|2)  =F{ho)  =  J ,/«.M)*M)MrW(A)  dr  d^ 

=  J ~th{T)K{T)  dr  =  7. 

We  can  conclude  that  in  order  to  maximize  the  SNR,  7  must  be  the  largest  eigenvalue  of 
the  kernel  Rq(t,  A)¥(t,  A),  and  /i£(t)  is  the  corresponding  eigenfunction. 

Assuming  perfect  synchronization,  i.e.,  p(p)  =  %),  it  is  interesting  to  point  out  that 
for  the  case  of  a  single  nonfading  path,  i.e.,  Rq{t,  X)  =  R  (constant)  for  all  r  and  A, 

7^o(T)  =  /  Ra(r,X)s0{T)s*0{X)h0{X)  dX 

=so(r)  J  Rsl{X)hl(X)  dX 
=Cs0(t), 

where  C  €  C.  This  is  the  familiar  matched  filter. 


(3.22) 


(3.23) 
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On  the  other  hand,  for  the  case  of  nonfading  multiple  paths, 

7/i*(T)  =  ^{|&|2  J  Rs*0{ A  -  &))MA)  ^A}s0(t  -  4>i) 

i 

=  ^2CiS0{T  ~  fa), 


(3.24) 


which  is  the  RAKE  receiver  [21].  The  derivation  presented  above  subsumes  and  generalizes 
these  important  special  cases. 


3.3  Optimal  Signaling 

We  would  like  to  extend  our  discussion  now  to  the  study  of  how  to  characterize  both 
the  transmitted  signal  s{t)  and  the  detection  filter  h(t)  so  that  for  a  given  fading  channel 
with  autocorrelation  function  Ra( r,  A)  =  E(a(r)a*(\)),  the  received  SNR  at  the  receiver 

is  maximized. 

For  a  fading  channel,  c(r;t)  whose  statistical  properties  are  known,  when  the  signal 
s{t)  is  sent,  the  received  waveform  can  be  represented  as 

r(t)  =  c(t;  f)  *  s(t  -  r)  +  n{t),  (3-25) 

where  n(t )  is  an  AWGN  with  power  spectral  density  N0/ 2.  Consider  the  problem  of 
designing  a  filter  h{t)  matched  to  c(r;t)  *  s(t),  then  shaping  s{t)  to  obtain  the  maximum 
SNR  possible.  That  is,  find 

r 

smax(*)  =  arg  max  |  mgc 

/  |s(t)|2dt=l  /  |/i(i)|2dt=l 

3.4  Characterization  of  Optimal  Signaling  for  Single-Path  Fading 

We  derive  the  CP  first  for  single  path  fading  with  no  synchronization  error.  This  will 
bring  us  better  understanding  of  the  properties  of  this  optimization  problem  (3.26).  We 
would  like  to  characterize  the  optimal  s{t)  and  h{t)  so  that  for  a  given  Rq(t,  A),  the  SNR 
at  the  receiver  is  maximized.  We  impose  the  constraint  that  both  s(t)  and  h(t)  have  unit 

energy, 

G(s)  =  j  s{t)s*(t)dt  =  1  and  H(h)  =  J  h(t)h*{t)dt  =  1.  (3.27) 


|/[c(r;t)  *  s(t)]h(t)dt 
No 


(3.26) 
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Without  synchronization  uncertainty,  the  received  signal  at  the  output  of  a  detector  filter 


=  J[a{t)s(t)  +  n{t)]h{t)dt, 


(3.28) 


with  the  corresponding  signal  and  noise  components  being 


rs  =  /  a(t)s(t)h(t)  dt ,  and  rn=  n{t)h{t)  dt. 


(3.29) 


Maximizing  the  SNR  is  equivalent  to  maximizing  the  ratio  between  the  variance  of  rs 
and  rn.  Because  of  (3.27),  we  have  £(|rn|2)  =  JV0/2.  So,  what  we  really  need  to  do  is  to 
maximize  the  variance  of  rs,  which  can  be  represented  as  the  following  functional 


£(|rs|2)  =EiJ J  aiT)s(T )h(r)a* (X)s*(X)h* (X)drdX] 

=  IJ  E[a{T)a*{X)]s{T)h(T)s*(\)h*{\)dTd\ 

=  U  Ra(r,X)s(r)h(r)s*(X)h*(X)drdX 


(3.30) 


=  F{s,h). 

Using  the  method  of  Lagrange  multipliers  to  maximize  F(s,  h),  i.e.,  we  construct  the 


Lagrangian 


J(s,h )  =  F{s,h )  -7o G(s)  —  7ii7(/i), 


(3.31) 


and  set  the  gradient  of  J{s,h)  at  (s0,/i0)  to  be  zero,  where  s0,h0  e  C2  are  the  extreme 
points  in  the  domain  of  F.  This  gives 


■§^J{so,ho)  fcF(s0,ho)  dsG(so) 

=  -  7° 

JrJ(so,M 


rsH(h0) 

mH(h0) 


(3.32) 


which  can  be  rearranged  to 


■§^F(sq,  h0)l  \&G(so)]  \mH(ho) 


=  70  a 
7n-F(so,ho)  Iq} 


(3.33) 


We  have,  from  Appendix  A,  that 


^C7(so)  =  4  ^G(so)  =  0,  §;H(ho)  =  °>  and  =  K- 


(3.34) 
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Setting  s{t)  =  x (t)  +jy(t)  and  h(t)  =  u(t)+jv(t),  we  obtain 


d  \  l{dF{s0,ho)  . dF{so,h0 ) 

T/(*o,  M  =  j  I  -  J  ^  - 


(3.35) 


-F(soM)  =  i{— 


1  fdF{s0,h0 )  .dFjsoM) 


Let  e  G  R,  F  €  £2,  we  obtain  the  following  derivative  results. 


(3.36) 


F(((x 0  +  eV)  +  jyo),ho)  -  F(s0,h0) 

=  JJ  Ra(r,  A)((x0  +  eV)  +  jyo)(r)h0{T){(x0  +  eV)  -  jyo){X)h*0(X)  dr  dX  -  F(s0,h0) 

«  J j  Rol{t,  X)[€V(t)sq(X)  +  eV(X)so{T)]ho(T)hl(X)  dr  dX 

=2e  f  f  Ra(T,X)ho(r)h*0{X)Re{s0{X))  dX  V(r)  dr, 

J  »/  /O  o»T\ 


(3.37) 


which  implies 


dF(s0,h0)  =  2  r  Ha(T)  A) /i0  (r) /iq ( A) Re (s0 ( A) )  dA.  (3.38) 

9  a:  J 


F(x0+j(yo  +  eV)),ho)  -  F(so,h0) 


=  J J  Ra(T,X)(xo  +  j{yo  +  eV))(T)ho(T)(xo  -  j{yo  +  eV)){^)ho(^)  dx  dX  F(so,ho) 

~JJ  Ra{T,X)[j{eV(T)s*0{X)  -  eV(X)so{T))]ho{T)h0{X)  dr  dX 
=2 e  f  f  Ra{T,X)h0{T)h*0{X)lm{s0{X))  dAF(r)dr, 

J  J  ( o  on\ 


(3.39) 


which  implies 


d  F(sq,  ho) 


(r)  =  2  /  Rq(t,  A)/io(r)h5(A)Im(s0(A))  dA. 


(3.40) 
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F(sq,  {(tto  +  eV)  +  jv o))  -  F(a0,  M 

=  U  Ra{r ,  A) s0 (t)  ((mo  +  +  ivo)(r)a5(A)((«o  +  eV)  -  jv 0)(A)  dr  dA  -  F(s0, 4) 

&  J  J  Ra{r,  X)[eV(T)h*0(\)  +  eV{\)h0{T)]s0{T)s*0{\) 

=2e  J  J  i?a(r>  A)so(r)s5(A)Re(/io(A))  dA  V(t)  dr, 


dT  dA 


which  implies 


(3.41) 


dF(so,ho)  =2fRa{Tt  X)s0(t)sU  A)Re(MA))  <»•  (3-42) 

aw  J 


F(sq,  (wo  +  j(vo  +  ^F)))  —  F(soi  ^o) 

=  J  J  Ra{r ,  A)(«0  +  j  (w0  +  eV)){r)ho  {t)(uq  -  j(v  0  +  eF))(A)/io(A)  dT  dA  -  F(s0,  h0) 
njj  iJa(r,A)[j(eV(T)/»5(A)-e7(A)fto(T))]«o(T)ao(A) 

=2e  J J  Ra(T,A)so(r)sS(A)Im(/i0(A))  dAF(r)dr, 


dr  dA 


so 


9F(s0,M(r)  =  2  f  Rq(t,  A)s0(r)so(A)Im(/io(A))  dA. 

aw  J 

Combining  the  results  from  (3.34), (3.38),  (3.40), (3.42),  and  (3.44),  we  obtain 


_a 

9s 


and 


■^rF(so,ho)  =  so(t)  f  Rq(t,  A)s$(A)h$(A)  dA. 
oh  J 


So,  (3.33)  implies 


'L(t)Mt)' 

70^0 (T) 

_L(t)s0(t)_ 

_7i^o(r)_ 

From  (3.47),  we  have 


L2(t)/i0(t)s0(t)  =  7o7i^o(t)so(t)> 


(3.43) 

(3.44) 


F(so,ho )  =  ho(r)  J  Ra(r,  A)sJ(A)/ig(A)  dA  —  ho(T)L(r),  (3.45) 


(3.46) 


(3.47) 


(3.48) 
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(3.49) 

(3.50) 

(3.51) 

(3.52) 


That  means  at  the  extrema  of  F{s,  h),  s(t)  and  h{t)  must  be  similar  in  magnitude. 
In  addition,  We  have  the  following  three  observations: 

1.  From  (3.47),  we  can  observe  that  L(r),  s0(r),  and  /i0(r)  have  the  same  zeros. 


That  means  the  SNR  of  the  CP  is  determined  by  the  square  root  of  the  product  of 
the  Lagrange  multipliers. 


3.  Let  {<pk}  be  the  set  of  eigenfunctions  of  the  symmetric  kernel  Ra  with  corresponding 
eigenvalues  {«&}, 

Kfc^fcO")  =  J  Ra{r,  A)^jk(A)  dX.  (3.54) 

Assume  that  the  eigenvalues  are  distinct  for  different  k.  Then  from  Tricomi  [115, 
Chapter  3], 

ih^k)  =  J  4>h{r)<l)k(T)dT  =  0  if  h^k. 

and,  span({^fc|)  =  C2{a,b),  with  (a,  b)  being  the  support  of  Ra. 


(3.55) 
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Now,  let  g( A)  =  So(A)ho(A)  —  Y^iak(t)k(^)i  then 
L(t)  =  j  Ra{r,  X)g(X)dX 

=  J  Ra{r ,  A)]Tafc0fc(A)dA  =  ^  akKk4>k{T)  • 


(3.56) 


Because  4>k  are,  orthogonal  and  nk  are  distinct,  we  should  be  able  to  characterize 
ofc  so  that  \L{t)\  =  constant.  That  means  we  are  able  to  find  /(A),  and  in  turn  the 
product  s*0{X)h*0{X).  Now,  because  of  (3.52),  we  may  determine  s0(A)  and  h0{ A)  up 


to  complex  conjugate  pairs. 
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Chapter  4 

Iterative  Algorithms  for  Finding  Solutions 

The  process  of  finding  CPs  s{t)  and  h{t)  directly  may  be  prohibitively  difficult.  In 
this  chapter,  we  present  iterative  algorithms  for  finding  these  CPs.  Prom  the  fact  that  the 
received  signal  component  is 

r  -  j  a(t)s(t)h(t)  dt,  (4-l) 

which  is  symmetric  in  s(t)  and  h(t),  if  we  are  able  to  characterize  h(t)  that  maximizes  the 
SNR  for  a  given  s(t),  we  are  also  able  to  characterize  s(t)  that  maximizes  the  SNR  for  a 

given  h(t). 


4.1  Iterative  Algorithms 

Assume  that  the  detection  filter  is  h0(t),  when  s(t)  is  transmitted.  Then  the  received 
signal  at  the  output  of  ho(t)  is 

ai(t)s(t  —  <pi)  +  n(t)]ho(t  +  p)  dt,  (4.2) 

and  the  signal  and  noise  components  are 

rs  =  J  -  <f>i)ha {t  +  dt,  (4.3) 


and 

rn  =  J  n{t)h0(t  +  n)  dt.  (4-4) 

Again,  we  have  E{\rn\^)  =  Nq/2.  The  variance  for  the  signal  part  is 


£;(|r5|2)  =E[J J (Y^&ai(T)s(T  ~  &))Mr  +  ^)(^^jaj(^)s*(^  +  v)  dr  dX] 

=  EE// ^jE[ai{T)a*{X)h0{r  +  +  h)}s{t  -  (j)i)s*{\  -  <pj)  dr  dX. 


(4.5) 


(4.8) 
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As  a(t)  and  so(i)  are  independent,  (4.5)  becomes 

j5(|rs|2)  =  ^  ^2  J J  £i€jE[ai(T)aj(X)]E[ho(T  +  n)hl(X  +  y,)]s(r  -  <f>i)s  (A  -  (f>j)  dr  dX. 

*  j  (4-6) 

Then,  the  assumption  of  independent  paths  gives 
E(\rs\2)  =  J J  \Ci\2E[ai(T)al(X)]E[ho(T+fi)hQ(X+fi)]s(T  —  (j>i)s*(X  —  <j>i)  dr  dX.  (4.7) 

Defining 

T  (r,  A)  =E[h0{T  +  li)h*0{X  +  n)} 

=  J  ho(r  +  h)Iiq(X  +  [i)p(n)  d/j,. 

Then  (4.6)  can  be  rewritten  as 

E(\rs |2)=X]  J j  |^|2^a(r,A)T(r,  A)s(r  -  &)s*(A  -  fa)  dr  dX 

—  J  J  RaiT  ~  fai  ^  ~  ^i){^  \(i\2~£(T  +  (f>i,X  +  (j>i)}s(T)s*(X)  dr  dX. 

Defining 

0(r,  A)  =  |&|2T(t  +  A  + 

i 

Again,  we  have  a  functional  expression  of  the  variance  for  the  signal  part. 

£(|rs|2)  =  / J  Ra(T,X)Q{T,X)s(T)s*{X)  dr  dX 


(4.9) 


(4.10) 


(4.11) 


=  F(s). 


We  use  the  method  of  Lagrange  multipliers  to  derive  s{t)  that  maximizes  F(s)  under  the 
constraint  of  G(s)  =  J  s2(t)  dt  =  1  by  setting 


dF(so)  _  dG(so) 
ds  ^  ds  ’ 


(4.12) 


where  so  €  £2  is  the  extreme  point  in  the  domain  of  F. 

For  e  6  K,  V  €  £2  with  s{t)  =  x{t )  +  jy{t),  we  have 
F((a;o  +  eV)  +  jyo)  -  F(so) 

=  J J  Ra{r,X)@{T,  X){(x0  +  eV)  +  jy){T){{xo  +  eV)-jy)(X)dTdX-F(so) 
kJJ  Ra(T,X)Q(T,X){eV(T)s*0(X)  +  eV(X)so(r)}  dr  dX 

=2e  J  J Ra(r ,  X )©(r,  A)Re(s0(A))  dA  F(r)dr 


(4.13) 
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both  Ra  and  0  are  symmetric  in  r  and  A,  and 
F(x o  +  j{yo  +  eV))  -  F(s0) 

—  j  j  Rol{t,X)Q(t,  X)(xq  +j{y  +  eV)){r)(xo  -  j(y  +  eV))(X)dTd\  —  F(sq) 

kJJ  Ra(T,X)e(T,X){j[eV(r)s*0(X)  -  eV(X)s0(r)}}  dT  dX 

=2e  J  J  i?a(r,A)0(r,A)Im(so(A))  dX  V{r)dT. 

That  is, 

dF(8°l(T)  =  j Ra(Ti  A)0(r,  A)sq(A)  dX. 


ds 


Equation  (4.15)  and  (4.12)  imply 

7«o(T)  =  J #a(r>  ^)©(T)  A)so(A)  dX. 

Additionally, 

£(|rs|2)  =F{s0)  =  j  j  i?Q(r,A)0(r,  A)s0(r)s5(A)  dr  dX 

=  J 7so(r)so(r)  dT 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


=7- 

Thus  to  maximize  the  SNR,  we  set  7  to  be  the  largest  eigenvalue  of  the  kernel  Ra(r ,  A)0(r,  A) 
and  set  Sg(f)  to  be  the  corresponding  eigenfunction. 

By  combining  the  results  from  (3.22)  and  (4.17),  we  obtain  an  algorithm  to  find  the 
CP  for  a  multipath  fading  channel.  That  is,  we  iteratively  fix  the  signal  function  s(t)  and 
the  detector  function  h(t),  then  evaluate  for  the  corresponding  optimal  filter  function  h0{t) 
and  signal  function  so(t).  This  is  summarized  in  the  following  algorithm. 


Algorithm  1  1.  Randomly  select  an  initial  function  so{t). 


2.  Set  ho{t)  so  that 

70 K(t)  =  J  Ra(r ,  A)tf  (r,  X)h*0(X)  dX,  (4.18) 

where  70  is  the  largest  eigenvalue  in  magnitude  of  the  kernel  i?a(r,  A)^(r,  A),  with 
■$(t,X)  =  ^  lCi|2^[so(r  -  f>i)s 5(A  -  f>i)]).  Then  h^ft)  is  the  eigenfunction  corre¬ 
sponding  to  70  • 
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3.  Set  si(t)  so  that  for  the  h0{t)  in  the  previous  step, 

71*0(1-)  =  / i?a(r,A)0(r,A)S5(A)  dX,  (4.19) 

where  71  is  the  largest  eigenvalue  of  the  kernel  i?Q(r,  A)0(r,  A),  with  0(r,  A)  = 

Si  |&|2-E![MT  +  <Pi)hoi^  +  </>i) ])•  Then  so(f)  is  the  corresponding  eigenfunction. 

4.  Repeat  step  2  to  step  4  a  sufficient  number  of  times. 

When  we  try  this  algorithm  with  discrete  samples  of  Ra,  s(t),  and  h(t)  with  the 
assumption  of  perfect  synchronization  using  Matlab,  we  find  that  |s(t)|2  and  \h(t)\2  both 
approach  5(f).  This  actually  makes  sense,  because  when  there  is  uncertainty  in  the  channel, 
one  would  like  to  transmit  information  through  that  channel  in  as  short  a  time  as  possible 
to  reduce  the  uncertainty  in  the  received  signal. 

Analogous  to  the  characterization  of  the  CP  for  multiplicative  channels  which  is  in 
the  time  domain,  we  may  derive  the  same  arguments  for  frequency  selective  channels  in 
the  frequency  domain.  This  amounts  to  having  a  multiplicative  factor  in  the  frequency 
domain.  In  this  case,  we  will  end  up  with  a  CP  as  5  functions  in  the  frequency  domain. 
This  agrees  with  the  observation  as  pointed  out  by  Biglieri  [22,  p.2636]  that  “...  some  recent 
insights  into  the  peaky  nature  of  capacity-achieving  signaling  in  the  realm  of  broadband 
time-varying  channels....”  In  addition,  a  well  known  capacity-achieving  signaling  method 
that  attracted  much  attention  is  frequency  shift  keying  (FSK),  which  is  “peaky”  in  the 
frequency  domain. 

4.2  Convergence 

The  precise  conditions  for  which  Algorithm  1  converges  are  unknown.  From  our  expe¬ 
rience  in  implementing  this  algorithm  with  discrete  approximations,  it  converges  most  of 
the  time,  at  least  in  the  cases  when  “exact”  eigendecompositions  are  used.  On  the  other 
hand,  proving  the  monotonic  increase  of  the  largest  eigenvalue  through  consecutive  itera¬ 
tion  of  Algorithm  1  for  the  case  of  single  path  fading  with  no  synchronization  uncertainty 
is  relatively  easy. 
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Theorem  1  The  sequence  of  the  largest  eigenvalue  generated  through  iterations  of  Algo¬ 
rithm  1  is  monotonically  increasing  for  single  path  fading  channels  with  no  synchronization 
uncertainty. 


Proof: 

For  single  path  fading  with  no  synchronization  uncertainty,  (4.18)  becomes 

70 K(t)  =  I Ra(T,X)so(r)s*0(X)K(X)  dX.  (4.20) 

Also,  (4.19)  becomes 

7i^(r)  =  J Ra(r,X)ho(r)h*0(X)s*0(X)  dX.  (4.21) 

Considering  that  (4.20)  and  (4.21)  are  repeated  consecutively  by  iterative  Algorithm  1,  we 
may  consolidate  these  two  expressions  into  one,  i.e., 


%3i(r)  =  J  dX,  (4.22) 


where  7fc  is  the  largest  eigenvalue  for  the  kth.  iteration,  and  gk  replaces  s  and  h  alternatively. 
Now,  we  have  from  (4.22)  and  (3.27)  that 


7 k=  J  9k{r)g*k(T)  dr 

=  JJ  Ra(T,X)gk-i{T)g*k_l{X)gk{T)g*k{X)  dXdr 
>JJ Ra(r,X)gk-i(T)gk_1{X)gk-2(T)gk_2{X)  dXdr 
=7Jfc— l- 


(4.23) 


This  completes  the  proof. 

We  are  interested  in  finding  CPs  that  maximize  the  SNR,  and  the  SNR  is  related  to 
the  largest  eigenvalue.  So,  the  variation  of  the  largest  eigenvalue  through  the  iteration  can 
serve  as  an  indicator  of  the  progress  of  the  algorithm. 
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4.3  Introduction  of  Frequency  Constraints 

From  section  4.1,  we  found  that  CPs  of  single  path  fading  channels  with  no  syn¬ 
chronization  uncertainty  are  <5  distributions.  However,  8  distributions  are  not  practical 
signaling  waveforms.  One  possible  approach  for  designing  practical  CPs  that  are  “opti¬ 
mal”  yet  producible  is  to  introduce  frequency  restrictions  to  s(t )  and/or  h(t).  In  addition 
to  requiring  s(t)  and  h(t)  to  be  unit  powered,  we  impose  restrictions  on  the  frequency 
spread  of  s(t)  and/or  h(t).  Let  S(f)  =  f(s(t))  and  H(f)  =  F(h(t))  denote  the  Fourier 
transforms.  We  introduce  the  functional  W4(s)  and  Wh(h)  as  measures  of  the  energy  of 
s(t)  and  h(t),  respectively  in  the  frequency  range  ( -f,f)Hz .  Specifically,  we  let 

w.W-  J'sws'Wdv 

=  /“[n(^)S(„)]s*M  dr, 

/oo 

[s(t)  *  sinc(27r/<)]s*(t)  dt  (4.24) 

-OO 

/oo  r  oo 

s*(t)  /  s(r)  sinc(27r/(t  -  r))  dr  dt 

-oo  J  — OO 

/oo  roo 

/  sinc(27r/(r  —  A))s(r)s*(A)  dr  dX, 

-OO  j  —  00 


and  similarly 


Wh(h)=J_  H(f)H*(f)  df 

/oo  roo 

/  sinc(27r/(r  —  X))h(r)h* (A)  dr  d\. 
-oo  J —oo 


(4.25) 


We  reconstruct  the  Lagrangian  in  (3.14)  and  (4.12) 


(1-/3) 


mF(so) 

MF(h  o)J 


+p 


£Ws(so)  +  £sWh(ho) 
[fhWs(so)  +  fhWh(h0)\ 


=  To 


&G{»  o) 


+71 


t,mo) 


.  (4.26) 


where  f3  G  [0, 1]  is  a  weighing  factor  between  the  SNR  and  frequency  containment.  As 


d 


Wa(s)  is  independent  of  h,  we  have  —  Ws(s0)  =  0.  Again,  let  s(t)  =  x(t)  +  jy{t),  then 
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evaluate 

Ws({x0+eV)+jy0)-Ws{s0) 

=2/  J  J  sinc(27r/(r  -  A))((x0  +  cV)  +  jyo){T)i(xo  +  eV)  -  jy){ A)  dA  dr  -  Ws(s0) 
«2/e  J  J sinc(27r/(r  —  A))[so(r)V(A)  +  So(A)V(r))  ^ 

=4 :/e  J  J  sinc(27r/(r  -  A))Re(so(A))F(r)  dA  dr, 


(4.27) 


and 


wa(x0  +  i(yo  +  ^))  -w,(«0) 

=2/  J  J  sinc(27r/(r  -  A))(aj0  +j(yo  +  eF))(r)(x0  -  j{y  +  eV))(A)  dA  dr  -  Ws(s0) 
«2/e  J  J  smc{2Tvf{T  -  \))[j(s*0(\)V(T)  -  s0{t)V{\))]  dA  dr 
=4/e  J  J  sinc(27r/(r  —  A))Im(so(A))V(r)  dAdr. 


We  have 


5^ 

5s 


Ws{s0)  =  2 f  J sinc(27r/(r  -  A))aJ(A)  dA  =  JT(t). 


Similarly,  we  have  ^W/^/io)  =  0,  and 


-zrWh{ho)  -2 f  [  sinc(27r/(r  -  A))/iq(A)  dA  =  M(t). 

oh  J 

As  a  result,  (4.26)  becomes 


(4.28) 

(4.29) 

(4.30) 


"(1-/9)/ i?a(r,  A)0(r ,  A)sq(A)  dA  +  /3*T(t)' 

7oso(T) 

(1  -  /3)  /i?a(r,  A)tf  (r,  A)/»J(A)  dA  +  /3M(r)_ 

7i^o(r). 

(4.31) 


If  we  use  the  iterative  approach,  we  have  the  following  algorithm: 
Algorithm  2  1.  Randomly  select  an  initial  function  so{t). 


2.  Set  ho{t)  so  that 

71ho(r)  =  J[{1  ~  P)Ra{r,  A )tt(r,  A)  +  2/3/  sinc(27r/(r  -  \))]h*Q{\)  dA,  (4.32) 

7!  is  the  largest  eigenvalue  of  the  kernel  [(1  -  P)Ra(r ,  A)1I,(t ,  A)  +  2/3 /  sinc(27r/ (r  - 
A))],  while  /ij)(t)  is  the  corresponding  eigenfunction. 
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3.  Set  so(t)  so  that  for  ha(t)  from  the  last  step, 

7oSo(t)  =  y  [(1  —  /3)Ra{r,  A)@(r,  A)  +  20 f  sinc(27 r/(r  -  A))]sg(A)  dX,  (4.33) 

7o  is  the  largest  eigenvalue  in  of  the  kernel  [(l-/3)i?Q(r,  A)0(r,  A)+2y0/sinc(27r/(r- 

A))]  while  Sg(<)  is  the  corresponding  eigenfunction. 

4.  Repeat  steps  2,  3  and  4  a  sufficient  number  of  times. 

Now,  Algorithm  1  is  a  special  case  of  Algorithm  2. 

4.4  Signal  Sets 

To  communicate  over  a  fast  fading  channel  whose  amplitude  may  change  sign  in¬ 
stantly,  amplitude  and  phase  modulation  are  inappropriate,  and  something  akin  to  orthog¬ 
onal  signaling  is  necessary.  We  must  therefore  find  at  least  two  sets  of  signal/filter  pairs 
(so(t),ho(t))  and  (si(t),  hi(t)),  such  that  under  the  channel  influence  aft),  Si(t )  and  hj(t) 
have  a  minimal  correlation  for  i  ^  j.  This  minimal  correlation  criterion  may  be  fulfilled  by 
a  direction  application  of  the  CPs  generated  with  Algorithm  2  in  pulse-position  modulation 
(PPM).  Another  more  general  approach  requires  an  extension  to  Algorithm  2. 

From  Algorithm  2,  we  know  that  the  minimum  correlation  between  hj(t)  and  Si(t) 
can  be  obtained  by  setting  hj(t)  to  be  the  eigenfunction  corresponding  to  the  smallest 
eigenvalue  while  we  generate  hi(t)  for  given  Si(t),  at  the  same  time  setting  Sj(t)  to  be  the 
eigenfunction  corresponding  to  the  smallest  eigenvalue  while  we  generate  st(t)  for  given 
hi(t).  However,  the  kernel  must  be  changed  to  enforce  frequency  containment.  So,  while 
finding  the  eigenfunctions  corresponding  to  the  smallest  eigenvalues,  we  use  the  kernel  of 
the  inverse  operator  (![•],  Appendix  B)  of  (4.29)  and  (4.30)  in  place  of  2/sinc(27r/(r  -  A)) 
in  (4.31)  for  frequency  containment. 

We  suggest  the  following  procedure:  Given  a  transmitted  signal  soft),  let  the  cor¬ 
responding  filter  ho{t)  be  the  eigenfunction  having  largest  eigenvalue  of  the  kernel  (1  - 
0)Ra(T,  A)’4r(r,  A)  +  2/?/sinc(27r/(r  -  A))  as  defined  before.  At  the  same  time,  set  hi(t) 
to  be  the  eigenfunction  with  smallest  eigenvalue  of  the  kernel  (1  —  /3)i?Q(r,  A)^f(r,  A)  + 
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/9J[2/sinc(27r/(T  -  A))].  Similarly,  given  a  detection  filter  hQ(t),  we  let  the  corresponding 
signal  so(t)  be  the  eigenfunction  having  largest  eigenvalue  of  kernel  (l-/3)f?Q(r,  A)$(r,  A)  + 
2/3/sinc(27r/(r  -  A)).  At  the  same  time,  set  s\(t)  to  be  the  eigenfunction  with  small¬ 
est  eigenvalue  of  the  kernel  (1  —  (3)Ra(T,  A)$(r,  A)  +  sinc(27r/(r  —  A))].  Use  the 

same  method  to  obtain  si(f),  hi(t),  so(t),  and  ho(t).  This  procedure  will  be  carried 
out  repeatedly:  for  each  new  iteration,  set  so(t)  =  so{t)  +  so(t),  si(f)  =  s\(t)  +  si(t), 
hQ(t)  =  ho(t)  +  ho(t),  and  hi(t)  =  hi(t)  +  h\(t),  then  normalize  these  functions  to  unit 
energy.  The  procedure  is  summarized  in  Fig.  4.1. 

- - ^  s0(t)  = 

^So(t)+So(t) 

/  M  «l(t)  = 

h(t)  +  h(t) 

ng  signal/filter  pairs. 

In  the  above  procedure,  so(^)  and  ho(t)  form  a  CP,  while  si(i)  and  hi(t)  form  another 
CP.  Further,  the  correlation  between  so(t)  and  hi(t),  and  that  between  si(t)  and  ho(t)  are 
small.  We  call  CPs  with  minimal  correlation  as  conjugate-consonant  pairs  (CCPs). 

Due  to  the  averaging  in  the  generation  process  of  the  signal/filter  pairs,  the  CCPs  vary 
a  lot,  both  with  the  waveforms  and  the  SNR  in  terms  of  the  largest  eigenvalue.  However, 
during  the  process,  we  come  across  many  “good”  pairs  which  have  high  SNR. 

Algorithm  3  1.  Randomly  select  initial  functions  so(t)  and  si(t). 

2.  Set  hoo{t)  so  that 

Too Kq(t)  =  J{(  1  -  P)Ra(r,  A)*o(t,  A)  +  2 /3f  sinc(27r/(r  -  A))]/&(A)  dX,  (4.34) 

7oo  is  the  largest  eigenvalue  of  the  kernel  [(1  —  /3)i?a(r,  A)^ro(/r)  A)  +  2/3/sinc(27r/(T- 
A))],  while  hoo(t)  is  the  corresponding  eigenfunction. 


so(t) 


^  h0(t)  = 

^  ho(t)  +  ho{t) 


•lit) 


hi{t)  = 
h{t)  +  hi(t) 


Fig.  4.1:  Procedure  for  generati 
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3.  Set  hoi  so  that 

foi*»oi(T)  =  J [(1  -/3)i?a(r,A)^0(r,A)  + /3I[2/sinc(27r/(r  -  A))]]^X(A)  dX,  (4.35) 

£n  is  the  smallest  eigenvalue  of  the  kernel  [(1  -/3)i?a(r ,  X)^>q(t,  A)+/3T[2/  sinc(27r/(r- 
A))]],  while  /iQi(f)  is  the  corresponding  eigenfunction. 

4-  Set  hn(t)  so  that 

7n hh(T)  =  fii1-  P)rc*(t,  A)$i(t,A)  +  2j8/sinc(27r/(r  -  A))]/iX0(A)  d\,  (4.36) 

7n  is  the  largest  eigenvalue  of  the  kernel  [(1  —  P)Ra(r ,  A)'5i(r ,  A)  +  2/3/  sinc(27r/(r  — 
A))],  while  h*n(t)  is  the  corresponding  eigenfunction. 

5.  hio  so  that 

£i0/4(t)  =  J [(1  -  /3)Ra{T,  A)^i(r,  A)  +  /31[2f  smc(2n f{r  -  A))]]^0(A)  dX,  (4.37) 

£10  is  the  smallest  eigenvalue  of  the  kernel  [(l—(3)Ra(T ,  A)'^i(r,  A)+/3X[2/  sinc(27r/(r— 
A))]],  while  h*0(t)  is  the  corresponding  eigenfunction. 

6.  Find  the  average,  i.e., 


ho(t)  =  hoo(t)  +  hio(t), 

(4.38) 

hi(t)  =  /ioi(i)  +  hn(t). 

(4.39) 

7.  Normalize  ho{t)  and  h\(t). 

8.  Set  soo{t)  so  that 

looser)  =  j [( 1  -  /3)i?a(r,  A)0o(r,  A)  +2y8/sinc(27r/(r  -  A))K0(A)  dX,  (4.40) 

700  is  the  largest  eigenvalue  of  the  kernel  [(1  —  /3)J?a(r,  A)©o(r,  A) +  2/3/  sinc(27r/(r  — 
A))],  while  Sqo(^)  2S  corresponding  eigenfunction. 
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9.  Set  soi  so  that 

CoisSi(T)  =  JK1  -  j8)i2a(r,  A)0o(r,  A)  + /3I[2/sinc(27r/(r  -  A))]]sqi(A)  dX,  (4.41) 

£u  is  the  smallest  eigenvalue  of  the  kernel  [(1— P)Ra(r,  A)@o(r,  \)+/31[2f  sinc(27r/(r— 
A))]],  while  SqiW  is  the  corresponding  eigenfunction. 

10.  Set  sn(i)  so  that 

(r)  =  J[(  1  ~  P)R«(t,  A)0i (t,  A)  +  2/3/  sinc(27r/(r  -  A))K0(A)  rfA,  (4.42) 

7ii  is  the  largest  eigenvalue  of  the  kernel  [(1  -/3)iia(r,  A)0i(r,  A)  +  2/3/sinc(27r/(r  — 
A))],  while  Sii(i)  is  the  corresponding  eigenfunction. 

11.  sio  so  that 

&o*;o(r)  =  /[(!"  A)0i(r,  A)  + /3T[2/sinc(27r/(r  -  A))K0(A)  d\,  (4.43) 

£io  is  i/ie  smallest  eigenvalue  of  the  kernel  [(1  -f3)Ra(T,  A)@i(r,  A)+/3I[2/sinc(27r/(r— 
A))]],  while  s*w(t)  is  the  corresponding  eigenfunction. 

12.  Find  the  average,  i.e., 


So(t)  =  SooOO  +  Slo(*), 

(4.44) 

Sl(t)  =  S0l{t)  +  Suit). 

(4.45) 

13.  Normalize  so(f)  and  Si(t). 

14-  Repeat  steps  2  to  13  a  sufficient  number  of  times.  Save  those  CCPs  with  high  SNR 


along  the  iteration. 
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Chapter  5 

The  Study  of  Probability  of  Error 

5.1  Introduction 

In  conjunction  with  the  algorithms  suggested  in  the  previous  chapter,  we  would  like 
to  develop  a  method  to  evaluate  the  performance  of  the  derived  CPs  in  terms  of  SNR 
versus  probability  of  error,  P(e).  So,  we  may  compare  the  performance  between  the  CPs 
generated  through  those  procedures  in  chapter  4  and  other  traditional  waveforms. 

Up  to  this  point,  we  have  not  imposed  any  specification  for  the  fading  process  a(t), 
so  the  principles  we  derived  can  be  applied  to  general  fading  channels.  However,  for 
the  following  derivation,  we  will  assume  the  fading  process  is  a  stationary  band-limited 
Gaussian  process  to  facilitate  our  development.  The  result  we  obtain  will  be  applicable  to 
Rayleigh  fading  channels,  which  are  the  most  common  fading  channel  models. 

5.2  Quadratic  Forms  of  Complex  Normal  Vectors 

Let  {vo, i>i,  •  •  •  , ujv-i}  be  complex  random  variables  with  real  and  imaginary  parts, 
vk  =  xk  +  jyki  which  are  normally  distributed  and  satisfy  the  following  relationships: 

E[{xi  -  xi){xj  -  Xj)]  -  E[(yi  -  yi)(yj  -  yj)]  (5.1) 

and 

E[{xi  -  Xi){yj  -  yj)}  =  -E[{xj  -  Xj)(yi  -  y;)]>  (5-2) 

where  Xk  —  E(xk)  and  yk  =  E{yk)-  Define  the  complex  normal  vector  v  as 

v0 

Vi 

.  ,  (5-3) 

UAT-l 


with  mean  E(v)  =  0,  and  covariance  M  =  E[wH],  A  quadratic  form  of  v  is  defined  as 

g  =  vhQv,  (5.4) 
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where  Q  is  any  Hermitian  symmetric  matrix. 

According  to  Turin  [116],  the  probability  density  function  of  g,  p{g),  has  a  character¬ 
istic  function 

*(‘)  =  nrrk’  <5-5> 

where  Qk  are  the  eigenvalues  of  MQ.  To  simplify  our  development,  we  assume  that  all  Qk 
are  distinct.1  Since  Qk  =  0  does  not  affect  $(£),  we  assume  that  all  the  eigenvalues  are 
nonzero. 


(5.7) 


The  probability  density  function  of  g  is 

1  f°° 

pig)  =  2 -  J  $(*)  exP(-i^)  dt, 

which  can  be  evaluated  by  using  the  residue  theorem, 

pig )  Jc  $(*)  exp i-i*g)  dz 

—j  ^{residues  enclosed  in  the  counterclockwise  contour}, 
where  C  is  a  contour  in  the  complex  plane  containing  the  real  axis.  The  integrand  in  (5.8) 
can  be  rearranged  as 

*(*)  exp i-jzg)  =  |  JJ  1  -  )zgn  }  exp(~Jzg^ 


(5.8) 


= ( n  -i  I  exp i-Jz9)- 

{^z  +  JQn  J 


(5.9) 


For  each  isolated  singular  point  —jei  \  we  can  represent  the  integrand  by  Laurent  series 
§(z)exp(-jzg) 

\ 

-1 


=  < 


N  .  _x 
JJ  JQn 


JQi 


N 


.  _i  t  ,  ny  exp(  jzg) 
Z+JQn  Z  +  JQi 


1  -  -  /  ■  •  -u  ,  ( z  +  3Qi X)2 


L  Z  +  jQ, 


—  {l  +  (z+jQi  )  + 


->]} 

(5.10) 


1When  there  are  repeated  eigenvalues,  the  characteristic  function  is 

•W-IIjpsgsr. 


(5.6) 


where  fin  is  the  number  of  replications  of  gn-  There  will  be  some  minor  modifications  to  our  derivation 
when  \in  7^  1  for  some  n.  To  simplify  our  discussion,  we  assume  there  are  no  duplicated  eigenvalues. 


Hence,  the  residues  of  $(2)  exp(— jgz)  are 


(at  1 

Res{$(z)  exp  {-jgz),  -jgj1}  =  jejj1  <  0t0  f  exP(“^r1^)  (5-n) 

n=l,njii  01  0n  ) 

A  contour  map  is  shown  in  Fig.  5.1.  In  this  figure,  there  are  two  possible  closed 
loops,  I\  + 12  and  I\  +  I3.  Which  contour  we  use  depends  on  the  sign  of  g:  e~9Z9  must  be 
bounded.  For  z  =  x  +  jy,  we  have  the  following  cases: 

1.  g  >  0:  e~9°z  =  e~"x  ■  e9y  implies  that  we  need  to  take  y  <  0  and  the  loop  C  =  I\  +I3. 
By  the  residue  theorem, 

N 

p{g)  =  X  ft1  S  II  o--  o'  r  exP(_^rl3)  (5-12) 

8i> 0  \n=l, n^i  01  0n  j 

2.  g  <  0:  e~99Z  =  e~99X  •  e9y  implies  that  we  need  to  take  y  >  0  and  C  =  I\  + 1^-  So, 


Fig.  5.1:  Contour  map. 
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5.3  Cumulative  Distribution  Function 


Theorem  2  Let  v  be  a  complex  normal  vector  with  mean  E(v)  =  0,  and  covariance 
M  —  E\yvH\.  Further  let  Q  be  any  Hermitian  matrix.  Then  the  quadratic  form  g  —  \HQ\ 
o/v  has  a  cumulative  distribution  function  of 

\  Ee,<o{n"1,„i«5A;}exp(-?erl)  if  (<  0 
5(0  =  \  1  -  £b>0  5*:}  -rf-ar1)  V  €20.  ■ 

w/iere  Qk  are  the  eigenvalues  of  MQ. 

Proof: 

The  computation  of  Fg{(),  the  cumulative  distribution  function  of  g,  falls  into  the  following 
two  cases,  following  (5.12)  and  (5.13). 

1.  (  <  0: 

Fg(0  =  f  p{g)  dg 

J  —  OO 


—  /  Ee.'1 


N 

,  n 

^  n=l,n^i 


Qi 


exp (-^  xg)  dg 


•  ft  Qn 


= -  E  «r‘ 


=  E 

£i<0 


JV 

n 


N 

n 

^ra=l,n^z 

ft 


ft 


•  ft  £?n 


f  exp  i-et  1g) 

J  — OO 


(5.15) 


dg 


exp(— *) 


^n=l,n^z 


.  ft  fti 


2.  £  >  0: 


rS(o = r  p(<?) 

J  — OO 


dg 


POO 

=1  -  y  p(fl)  dg 

roo  {  N 

-l~Jt  E «r‘ 

i-Eft-f  n 


n 


ft 


^n=l,n^z 


Qi  Qn 


>  exp(-^  dg 


(5.16) 


ft 


0t>O  ^n=l,n^z 

=i-e(  n 


■  ft  £n 
ft 


roo 

J  exp(-p“15)  dg 


Qi >0  I  n=l,n^z 


■  ft  £?n 


V  exp(-^Pj  4). 
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5.4  A  Partition  of  Unity 

The  following  results  do  not  relate  directly  to  our  development  of  CPs.  We  present 
them  here  as  they  are  very  interesting  and  also  very  useful  in  the  field  of  digital  com¬ 
munication.  When  £  =  0,  we  have  the  following  corollary  immediately  from  (5.15)  and 
(5.16). 


Corollary  1  Let  {01,  £>2)“*  ,£at}  be  distinct  real  numbers ,  then 

N  N 


e  n 


Qi 

Qi  ~  Qn 


=  1. 


(5.17) 


The  above  derivation  of  the  cumulative  distribution  function  proves  that  this  corollary  is 
true  for  real  sequences  of  {pi,  Q2,  ■  ■  •  ,  Qn}-  This  result  can  be  extended  to  include  distinct 
complex  sequences  [117]. 

The  above  partition  of  unity  can  also  be  expressed  by  the  following  beautiful  forms: 


1.  Let  Qi  =  igi,  then 


N  N  N  N 

e  n  =—e  n  rh-1 

i=l  n=l3n^i  71  2=1  n=l,n/i 

1  VV  UN~l  f  N  ~  1  ^  j^-l  —  1 


2.  When  Qi  =  \q\,  then 


N  N 


N  N 


e  n  ^=e  n  £-t-1 

2=1  n=l,n^i  2=1  n=l,n^i 


(5.18) 


(5.19) 


3.  For  =  x1  where  x  ^  {0, 1}, 


N  N 

e  n 

2=1  n=l,n^2 


(5.20) 
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4.  For  Qi  -  ix,  x  ±  0, 


N  N 


iY  1  *  ^ 

^  II  i  _  (n/i)x 


2=1  n=l,n^i 


(5.21) 


Other  than  being  pretty  in  expression,  the  above  partition  of  unity  can  also  be  used 
to  simplify  many  expressions  presented  in  the  literature,  for  example  [21,  p.802(14-5-27)] 
and  [92,  p.274(9.26)]. 


5.5  Quadratic  Receiver 

A  quadratic  receiver  is  a  receiver  that  employs  a  filter  bank  to  obtain  a  received  signal- 
plus-noise  vector,  then  determines  which  signal  is  transmitted  out  of  a  finite  number  of 
possible  signals  by  considering  the  signal-plus-noise  vector  as  a  whole.  Assume  the  channel 
is  C(t,  t )  with  AWGN  and  may  be  either  a  fading  or  a  nonfading  channel.  There  are  K 
possible  transmitted  signals  {s0{t),  si,  •  •  •  ,  SK-i(t)},  and  there  are  N  filters  in  the  receiver 
filter  bank  {h0{t) ,  h{t) ,  •  •  ■  ,hN-i(t)}.  Fig.  5.2  shows  a  quadratic  receiver.  Let  rk{t)  be 
the  received  signal  before  the  receiver  when  sk(t)  is  transmitted  through  channel  C(r,t ) 
with  AWGN  n(t).  Then 

rk(t)  =  c(t,  t )  *  sfc(t)  +  n(t),  (5.22) 


where  n(t )  has  a  power  spectral  density  of  iV0/2.  Define 
Tki  =  J rk{t)hi(t)  dt 

=  J  c(r,t)  *  sk(t)hi(t)  dt  +  J  n(t)hi(t)  dt 

=  J (f)sfc(t-  <f>j)hi(t)  dt  +  J n(t)hi(t)  dt, 


The  signal-plus-noise  vector  when  sk(t)  is  sent  is 

rk0 

Dfci 

rkN- l  J 


(5.23) 


rjfe  = 


(5.24) 
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Let  rjk  =  E(rk).  The  covariance 
defined  as 


matrix  of  the  sample  vector  when  sk{t)  is  transmitted  is 


Rk  =E[{  rk-rk){vk-rk)H] 

=E{ rferf  -  vkfkH  -  fkr%  +  rfcrf ) 

=E{ rfcrf )  -  E{ rk)rkH  -  rkE{r%)  +  ffeff 
=E{rkr%)  -fferf. 


(5.25) 


We  have 


rfc  = 


^ko 
?ki 

L  ffeN_i  J 


(5.26) 


where 


rk.  =  J  E[C(r,t )]  *  sk(t)hi(t)  dt  +  j  E[n(t)hi(t) ]  dt, 


(5.27) 


which  will  be  zero  for  fast  fading. 

To  simplify  our  discussion,  we  assume  that  K  =  2,  and  let 


Cfc  = 


n0rko  rkorkl 
[nAo 


(5.28) 


which  is  zero  for  fast  fading.  Also, 


E{  vkr”)  =E( 


n0r*ko  rko  rkl 
L  r*Ao  r^rki 


ri too  rfeoi 

No 

1 

r 

j  h0{t)hi{t)  dt 

_  rk io  ’"fell  . 

2 

/  h\{t)ho{t)  dt 

1 

(5.29) 
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where 


(Mr)sdr  -  *  -  rt)Wr)  *  / WK*  "  *»  +  ^ 

=B[/  /(E««xM>‘(t  -  *  -  '‘»<E«“i(XK(x  -  *  -  '‘))Wt)'‘5W  dT 
=EE  £i/  /  -  ♦i-rtW'W’ dT  dX1 

=  ^  J  j  Elai(T)a‘’jW}E[(i(jsUr  -  <tii  -  -  h  -  ^h’>^ho^  iT  dX 

=  E  /  /  -Rof1-"  >‘)'B[l€»|2«fc(T  -<#■<-  M)f t(A  -  -  J*)]ho<T>fe5(A)  <Jt  dA 

=  J  J  H„It,  X)  J^{|fil2E[*»(T  I»K0  -  "  f*)]}MT)Mx) *■  £iX 

=  [  [  JSa(r,A)^fc(r,  \)ho{T)h*0{\)  d,T  d\, 

J  J  (5.30) 


where  we  define 


Similarly, 


r,  A)  =  ^{|&|2£Mr  ~fa~  ^S*^X  ~  fa  ^ 


rfc01  =JJ  Ra{r,  A)^fc(r,  A)/i0(t)J»J(A)  dr  dA, 
rfcl0  =JJ  Ra{r ,  A)*fc(r,  A)/ii(r)h5(A)  dr  dA,  and 
rfcll  =JJ  Ra(r ,  A)«fc(r,  A)fci(r)M(A)  dr  dA. 

For  the  perfect  synchronization  situation, 

tffc(r,  A)  =  “  <&)sfc(A  -  &)• 

i 

Without  the  assumption  of  synchronization, 

*,(T,X)  =  E«?  /  «o(T  -  *  -  )*)*5(x  -  A  - 

i 

where  p(/x)  is  the  density  function  of  jx.  So, 

[  i  [  h0{t)hx(t)  dt  1 

f  7-i.nn  rnil  .  ■‘*0  d  - 


_  rfc0o  ’"fcoi  +  . 

Kfc  _  [  r*i°  rfcU  J  2  /  hi{t)h0{t)  dt 


-Ck. 


(5.31) 


(5.32) 
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Since  correlation  detection  is  a  linear  operation,  r  is  complex  Gaussian  with  density 
function  [118] 


P(r\sk(t))  =  exPHr  -  n)HRkl(r  -  **)]■ 


(5.33) 


Assume  that  both  signals  occur  with  equal  prior  probability.  The  maximum  likelihood 
decision  rule  we  use  is 


p(r\s0(t))  -  p{r\si{t))  <0 

1  50 

»^exp[-(r  -  ft)”  V(r  -  *»)]  -  «‘PHr  -  ^XR^r  ~  fl)1 1 ° 


?  |J*ol 


(5.34) 


=*(r  -  *i)*«[Hr  ~  fi)  -  (r  -  fo  W('  -  fo)  <  ln\R[\ 

Mr  -  u )HQ(r  -  u)  i  In  j^j  +  ff  R?f. i  -  ^  Rfy  ~  w  =  * 

so 

=>g  =  vhQv  <  £, 

Si 

where  Q  =  H?1  -  and  v  =  r  -  u.  Let 

Mo  =  S[(v-v)(v-v)H|so(t)]  and  Mi  =  J5[(v - v)(v - v)H|si(t)].  When  «o(*)  is  sent’  we 
havev|So(t)  =r  +  Q'1(^o1ro--Rr1ri)  and v|So(t)  =  f0  +  Q"W l*o -Ri^i)-  Therefore, 
(v  _  v)|so(t)  =  (r  -  f0)  which  implies  M0  =  Ro-  Similarly,  when  «i(t)  is  sent,  we  have 

Mi  =  Rv 

When  we  assume  fast  Ming,  n  =  0,  (5.15)  and  (5.16)  can  be  utilized  to  evaluate  the 
performance  of  different  signaling  waveforms.  For  {  >  0,  the  bit  error  probability  when 

S0(t)  is  sent  can  now  be  expressed  as 

K 
-00 


P(e|s0(t))  =  [  p(g\so{t))  dg 
7-00 

=Fg\so{t)(& 

=i-e(  n 


Qj{ 


di >0  ^n=l,n^i 


.  Cii  —  dn 


>  exp(-^ai  *), 


(5.35) 
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where  an  are  the  eigenvalues  of  RqQ.  The  bit  error  probability  when  «i (t)  is  sent, 

r  OO 

F(e|si(t))  =  /  p(fl|«i(*))  d9 

J  f 

=1  -  [  p(g\s\{t))  dg 

J  —  OO 

=1  —  Fg\si(t)(0 


=E 

6i>0 


N 

n 

n=l, 


.bi  —  bn 


^  exp(-^6j  *), 


where  &„  are  the  eigenvalues  of  R\Q-  Similarly,  for  C  <  0,  we  have 


N 


w.(«) = e  {  n  ^hpK°r) 

ai<0  ln=l,n^i 


and 


N 


?(ewt))  =  i-El  n  ^rxp(-e6rl)- 


bi<0  ln=l,n^i 


The  overall  bit  error  probability  is 


P(e)  =  ip(c|a0(t))  +  |P(e|»i(t))- 


This  expression  will  be  used  to  characterize  the  performance  of  the  waveforms 
using  the  iterative  algorithm. 


(5.36) 


(5.37) 


(5.38) 


(5.39) 

created 
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Chapter  6 

Implementation  of  the  Algorithms 

As  mentioned  in  section  1.3,  we  are  working  with  physically  realizable  signals  and 
filters  in  £2.  The  difficulty  in  finding  CPs  is  that  solutions  must  be  obtained  by  solving 
integral  equations.  These  depend  on  the  channel  statistics  properties.  As  is  frequently  done 
in  numeric  work  in  general  [119],  we  work  with  sampled  versions  of  these  functions.  Instead 
of  working  with  the  infinite  dimensional  functions  directly,  we  use  their  approximations  in 
finite-dimensional  spaces  to  be  able  to  implement  these  algorithms  from  chapter  4.  As  the 
sampling  interval  decreases,  a  solution  of  increasing  accuracy  is  obtained. 


6.1  Conversion  of  a  Linear  Integral  Equation  to  Its  Discrete  Equivalent 

A  linear  integral  equation  may  be  converted  to  a  discrete  equation  by  the  following 

approximation.  For  every  r, 

rb 

y(r)  =  /  R(t,  A)x(A)  dX 
J  a 


(6.1) 


JD 

=  lim  ^  s  .R(t,  mAA)x(7nAA) AA, 

AA->0  ^ 


m—A 


where  J,B£N  such  that  a  *  AAA  and  6  «  BAX.  So,  tor  AA  small  enough,  we  have 


B 


y(r)  «  ^  _R(r,7nAA)x(mAA)AA. 

m=A 

Now,  we  can  represent  y(r)  in  its  discrete  form. 

y(r)  =  lim  y(nAr), 

’  Ar-hO 

where  n  G  N  such  that  r  =  nAr.  Again,  when  At  small  enough,  we  have 

y(r)  wy(nAr). 


(6.2) 


(6.3) 


(6.4) 


Let  AT  =  min(Ar,  AA),  and  equating  (6.2)  and  (6.4),  we  obtain 


B 

y(nAT)  =  fl(™AT,mAT)x(mAT)AT. 

m—A 


(6.5) 
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Define  the  following  vectors  and  matrix  with  elements  being  the  values  of  y (nAT) ,  x(mAT), 
and  R{nAT,mAT),  respectively: 


r 

x((n  —  1)AT) 

y((n-l)AT) 

x(nAT) 

,  y  = 

y(nAT) 

x((n  +  1)AT) 

y((n  +  1)AT) 

_ 

i- 

and 


R  = 


R(mAT,nAT)  ••• 


(6.7) 


Then  (6.5)  is  equivalent  to 


y  =  RxAT. 


(6.8) 


Thus  the  integral  equation  is  approximated  by  a  matrix  equation.  We  will  apply  this 
conversion  to  the  pertinent  equations  from  chapter  4.  Equations  (4.32)  and  (4.33)  have 
the  following  form 


y(r)  =  f[Ra{r,  X)C{t,  A)  +  fcsinc(27r/(r  -  A))]y(A)  dX 

J  (6.9) 

=  J  Ra{r,X)C(T,X)y(X)  dX  +  J  fcsinc(27r/(r  -  A))y(A)  dX. 

In  (6.9),  we  have  the  discrete  equivalent  of  /  i?a(r,  A)C(r,  A)y(A)  dX  being  Py,  where 


•  •  •  Ra(mAT,  nAT)C{mAT,  nAT) 


=  RaoC.  (6.10) 


In  (6.10),  Ra  is  the  matrix  representations  of  Ra(r,  A),  C  is  the  matrix  representation  of 
C{t,  A),  and  o  represents  the  Hadamard  (element  by  element)  product.  So,  (6.9)  can  now 
be  represented  in  discrete  form  as 


y  =  [Rq  °Cf  fcW]y, 


(6.11) 


where  W  is  the  matrix  representation  of  2/ sinc(27r/(r  -  A)). 
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6.2  Discrete  Version  of  Algorithm  2  and  3 

Let  Ra,  $,0,  W,  and  S  be  the  matrix  representation  of  Ra{r,  A),  ^(r,  A),  0(r,  A), 

2 /  sinc(27r/(r-  A)),  and  X[2/  sinc(27r/(r  -  A))],  respectively.  Algorithm  2  can  be  presented 
in  sampled  form. 

Algorithm  4  1 ■  Randomly  select  an  initial  vector  So  of  length  N. 

2.  Set  ho  to  make 

7ih0  =  [(l-/3)RQo$  +  /rVV]h0,  (6-12) 

where  71  is  the  largest  eigenvalue  in  magnitude  of  the  kernel  [(1  /3)Ra'fr  -b  ^AV], 
while  ho  the  corresponding  eigenvector. 

3.  Set  s0  so  that  for  h0  from  the  last  step, 

7oSo  =  [(l-/3)Ra°0  +  £W]so,  (6.13) 

where  70  being  the  largest  eigenvalue  in  magnitude  of  the  kernel  [(1  -  /3)Ra©  +  (3W], 
and  §0  being  the  corresponding  eigenvector. 

4.  Repeat  steps  2,  3,  and  4  a  sufficient  number  of  times. 

Algorithm  3  can  be  presented  in  discrete  sampled  form. 

Algorithm  5  1.  Randomly  select  initial  functions  so  and  Si- 

2.  Set  hoo  so  that 

700 h00  =  [(1  -  j9)Ra  °  *0  +  /3W]h00,  (6.14) 

700  is  the  largest  eigenvalue  of  the  kernel  [(1  —  /3)Ra  0  ^0  +  /3W],  while  hoo  the 
corresponding  eigenfunction. 

3.  Set  hoi  so  that 

£oih0i  =  [(1  -  j3)Ra  o  *  0  +  iSSjhoi,  (6.15) 

£11  is  the  smallest  eigenvalue  of  the  kernel  [(1  —  /5)Ra  0^0  +  $S],  while  hoi  is  the 
corresponding  eigenfunction. 
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4 .  Set  hn  so  that 

Tiihn  =  [(1  -  £)Ra  o  *  i  +  /3W]h10,  (6-16) 

7n  is  the  largest  eigenvalue  of  the  kernel  [(1  -  /3)Ra  o  +  0W],  while  hn  is  the 
corresponding  eigenfunction. 

5.  Set  hio  so  that 

£iohio  =  [(1  -  jS)Ra  o  *  i  +  j8S]h10,  (6-17) 

£10  is  the  smallest  eigenvalue  of  the  kernel  [(1  -  /3)Ra  0  *i  +  whlle  ^io  IS  the 
corresponding  eigenfunction. 

6.  Find  the  average,  i.e., 

ho  =  hoo+hio,  (6.18) 

and 

hi  =  hoi+hn.  (6.19) 

7.  Normalize  ho  and  hi. 

8.  Set  Soo  so  that 

Toosoo  =  [(1  -  0)R a  °  ©o  +  i0W]soo,  (6-2°) 

7oo  is  the  largest  eigenvalue  of  the  kernel  [(1  -  /3)Ra  0  @o  +  /3W],  while  s00  is  the 
corresponding  eigenfunction. 

9.  Set  Soi  so  that 

CoiSoi  =  [(1  -  £)Ra  o  00  +  j8S]soi,  (6-21) 

£u  is  the  smallest  eigenvalue  of  the  kernel  [(1  -  /3)Ra  0  ©o  +  £S],  while  soi  is  the 
corresponding  eigenfunction. 

10.  Set  sn  so  that 

7uSn  =  [(1  -  /3)Ra  °  ©i  +  /3W]§io,  (6.22) 

7u  is  the  largest  eigenvalue  of  the  kernel  [(1  -  /3)Ra  o  ©i  +  /?W],  whlle  sn  %s  the 
corresponding  eigenfunction. 
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11.  Set  sio  so  that 

eiosio  =  [(1  -  0)Ra  °  ©1  +  /5S]sio,  (6-23) 

f10  is  the  smallest  eigenvalue  of  the  kernel  [(1  -  /3)Ra  0  0i  +  /3S],  while  sl0  is  the 
corresponding  eigenfunction. 

12.  Find  the  average,  i.e., 

so  =  soo  +  sio,  (6.24) 


and 


Si  =  s0i  +  sn. 


(6.25) 


13.  Normalize  So  and  Si. 

14.  Repeat  steps  2  to  13  a  sufficient  number  of  times.  Save  these  CCPs  with  a  large  SNR. 

As  Ra (r,  A)  is  an  autocorrelation  function,  while  2/  sinc(27r/(r- A))  and  l[2f  sinc(27r/ 
(j  _  /\))]  are  the  inverse  Fourier  transform  of  non-negative  functions,  all  these  kernels  are 
positive  semi-definite.  So,  Ra,  W,  and  S  are  required  to  be  positive  semi-definite,  too. 
Additionally,  since  they  are  Toeplitz  matrices,  the  requirement  of  positive  semi-definite 
implies  the  discrete-time  Fourier  transforms  of  the  sequences  of  their  elements  have  to  be 
non-negative  [120].  In  general,  a  sampled  sequence  is  not  necessary  non-negative.  This  is 
because  when  we  perform  discrete  sampling  from  a  continuous  function,  g,  and  present  it 
as  a  finite  dimensional  vector,  we  introduce  two  artifacts.  These  are  the  windowing  effect 
as  we  truncate  the  function  to  a  finite  length,  and  the  aliasing  effect  if  we  do  not  sample 
fast  enough.  Aliasing  can  be  controlled  fairly  well  since  most  physical  autocorrelation 
functions  are  low-pass.  When  we  perform  any  kind  of  windowing,  we  equivalently  convolve 
the  Fourier  transform  of  g  with  the  Fourier  transform  of  the  windowing  function.  By 
performing  this  convolution,  we  introduce  negative  values  to  the  Fourier  transform  of  the 
resulting  function.  The  effect  of  introducing  negative  values  to  a  positive  sequence  can 
only  be  reduced,  but  cannot  be  completely  eliminated,  if  we  use  a  windowing  function  that 


has  a  long  support. 
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The  sampled  sequence  can  be  made  to  have  a  positive  transform  again  by  running 
it  through  a  composite  mapping  algorithm  [121].  However,  this  usually  introduce  a  delta 
function  at  g(Q),  which  causes  a  deviation  of  the  sample  vector  from  the  original  function. 
Fig.  6.1  illustrates  this  deviation  and  the  new  peak  to  g( 0). 


6.3  Multipath  Fading  Channel  Model 

We  use  the  autocorrelation  function  of  the  land-mobile  model  (as  shown  in  Table  2.1) 
Rq(t)  =  j0(27rBdr)  to  represent  the  fast  fading  aspect,  and  the  European  Global  System 
for  Mobile  Communications  (GSM)  typical-urban  (TU)  and  hilly-terrain  (HT)  models  to 
represent  the  multipath  aspect  of  our  study  [96],  Table  6.1  presents  the  amplitude  and 
delay  profiles  of  the  simpler  6-path  TU  and  HT  multipath  channels. 


Table  6.1:  Delay  and  Amplitude  Profile  of  GSM  Channels. 


JiC  U  •  -L  • 

Path  # 

TU  model 

HT  model 

Delay  (/is) 

Power 

Delay  (/us) 

Power 

j 

0.000 

1.000 

0.000 

1.000 

2 

0.813 

0.669 

0.813 

0.251 

3 

1.626 

0.448 

1.626 

0.060 

4 

2.439 

0.300 

15.447 

2.258 

5 

3.252 

0.200 

16.260 

0.177 

6 

4.056 

0.134 

17.073 

0.122 
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6.4  Numerical  Methods 

The  most  prominent  challenge  we  face  in  computing  the  numerical  solution  to  (6.12), 
(6.13)  and  the  similar  equations  in  Algorithm  5,  is  the  size  of  the  matrices.  We  need  to  solve 
for  the  eigenvectors  corresponding  to  the  largest  eigenvalue  and  the  smallest  eigenvalue  of 
these  matrices.  A  large  matrix  size  exists  because  we  need  a  large  sample  size  to  take  off 
the  original  function  before  truncation.  For  example,  we  need  a  span  of  at  least  0  to  3 
seconds  for  J0(-)-  Only  then  the  sampled  sequence  resemble  the  original  sequence  loosely. 
This  is  shown  in  Fig.  6.2.  This  is  similar  to  Fig.  6.1,  which  is  shown  at  a  different  time 
scale.  We  can  observe  the  same  delta  function  at  time  0  but  with  a  smaller  magnitude. 
Also  shown  in  Fig.  6.2,  half  of  the  the  sample  sequence  is  set  to  zero  in  order  to  get  a 
proper  autocorrelation  matrix.  So,  the  actual  approximation  is  done  for  the  range  oft  from 
0  to  3  seconds.  Now,  some  simple  arithmetic  can  show  us  the  magnitude  of  the  sample 


Fig.  6.2:  Approximation  of  Jo{t)  with  a  positive  sequence. 

size  we  need.  From  Table  6.1,  we  can  see  that  the  path  delays  are  multiples  of  0.813/js.  It 
is  reasonable  for  us  to  use  one-tenth  of  the  unit  of  the  path  delays  as  our  sampling  time, 
i.e.,  ST  =  8.13  x  1CT8s,  where  ST  is  the  sampling  interval.  So,  We  have 

2irBd8Tn  =  3  =>Bdn  =  5872876.129,  (6-26) 

where  n  is  the  sample  number  corresponding  to  t  =  3  seconds,  Bd  is  the  Doppler  spread 
and  is  defined  as  Bd  =  j,  where  v  is  the  travel  speed,  and  A  is  the  wave-length  of  the 
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carrier  waveform.  Suppose  we  are  walking  at  a  speed  of  2  meters  per  second  while  talking 
on  a  cellular  phone  with  a  frequency  of  900MHz.  Then,  Bd  =  6Hz.  This  will  mean  we 
need  a  sample  size  of  978,813!  This  number  corresponds  to  a  978,813  x  978,813  matrix 
whose  eigenvalues  would  have  to  be  computed.  We  may  reduce  the  number  by  changing 
the  setting.  For  example,  if  we  are  traveling  in  a  fast  car  at  80  miles  per  hour  while  talking 
on  a  cellular  phone  with  a  frequency  of  1800MHz,  we  have  Bd  =  216Hz.  Furthermore, 
if  we  use  6T  =  1.626  x  10_7s,  we  need  a  sample  size  of  13,595.  Although  this  is  a  much 
more  implementable  figure,  we  still  need  to  improve  the  method  of  numerical  computation 
so  that  the  algorithm  may  be  implemented  on  today’s  computers.  Mainly,  we  want  to 
reduce  the  memory  requirement  and  eliminate  unnecessary  computations.  For  example,  the 
eigenvectors  corresponding  to  the  largest  eigenvalues  can  be  calculated  by  using  the  power 
method.  As  the  kernel  of  the  linear  equation  is  symmetric,  the  eigenvectors  corresponding 
to  the  smallest  eigenvalues  can  be  calculated  by  using  the  conjugate  gradient  method  [122]. 
So,  we  can  avoid  having  to  compute  all  the  eigenvectors. 

The  power  method  is  an  iterative  method  that  repeatedly  multiplies  an  initial  random 
vector  by  the  matrix  for  which  we  intend  to  find  the  eigenvector  corresponding  to  the 
largest  eigenvalue  in  magnitude  (which  may  be  either  positive  or  negative).  By  this  process, 
the  product  of  the  matrix-vector  multiplication  converges  to  the  eigenvector  with  largest 
eigenvalue  in  magnitude.  The  method  requires  that  there  be  a  single  eigenvalue  of  largest 
magnitude.  The  convergence  rate  of  the  power  method  is  highly  dependent  on  the  relative 
magnitudes  of  the  eigenvalues.  In  the  extreme  case  that  all  the  eigenvalues  have  the  same 
magnitude,  the  power  method  will  not  converge. 

The  conjugate  gradient  method  is  an  iterative  method  that  solves  an  equation  of  the 
form 

Ax  =  b,  (6.27) 

which  is  recommended  by  Demmel  as  the  algorithm  of  choice  for  systems  involving  sym¬ 
metric  positive  definite  matrices  [122].  What  this  method  does  is  actually  find  x  that 
minimizes  the  function 

/(x)  =  ^xAx  “  bx. 


(6.28) 
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When  the  above  function  is  minimized,  the  gradient 

V/(x)  =  Ax  -  b  (6.29) 

is  zero.  The  conjugate  gradient  method  itself  does  not  solve  the  eigenvalue  problem.  By 
iteratively  applying  the  conjugate  gradient  method  to  an  initial  random  vector,  we  have 
the  inverse  power  method.  The  process  converges  to  the  eigenvector  corresponding  to  the 
smallest  eigenvalue  in  magnitude.  The  convergence  rate  of  the  conjugate  gradient  method 
again  depends  on  the  relative  magnitude  of  the  eigenvalues  of  the  matrix.  The  closer  the 
small  eigenvalues  are  to  each  other,  the  slower  the  convergence. 

If  we  examine  the  expression  of  the  equations  in  Algorithms  4  and  5  more  closely,  they 
are  of  the  following  form 

7U  =  [/ 3R  o  (vv^)  +  (1  -  /3)Q] u,  (6.30) 

for  some  vector  v.  Since  R  and  Q  are  Toeplitz  Hermitian  matrices,  the  entire  computation 
can  be  done  with  vector  manipulations.  This  can  drastically  reduce  the  storage  used. 
Every  nxn  matrix  can  be  represented  using  an  n  x  1  vector. 

In  fact,  both  the  power  method  and  the  conjugate  gradient  method  only  require  being 
able  to  compute  matrix-vector  multiplications,  i.e.,  we  need  to  be  able  to  calculate 

[R  o  (vvH)  +  Q]x  =  y,  (6.31) 

assuming  that  f}  and  (1-/3)  are  absorbed  into  R  and  Q,  respectively.  Let  r  be  the 
coefficients  of  R,  and  q  the  coefficients  of  Q.  For  (6.31),  each  term  of  y  can  be  evaluated 
as 

Vi  =  Y^[r(j-i)viv*j  +  (6’32) 

3 

The  power  method  and  the  conjugate  gradient  method  in  cooperation  with  the  above 
method  of  matrix-vector  multiplication  are  used  to  evaluate  Algorithm  4.  A  result  is 
shown  in  Fig.  6.3  with  8192  =  213  sample  points.  Fig.  6.3,  like  many  other  plots  in 
chapter  7,  is  arranged  in  groups  of  four.  The  top  left-hand  plot  depicts  the  initial  random 
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signal  while  the  top  right-hand  plot  depicts  the  change  of  the  eigenvalue  through  the 
iterations.  The  bottom  left-hand  plot  depicts  the  final  s  while  the  bottom  right-hand  plot 
depicts  the  final  h.  There  are  two  curves  in  the  plot  of  eigenvalues.  One  of  the  curves 
is  obtained  from  the  iteration  of  the  s  part  while  the  other  curve  is  obtained  from  the 
iteration  of  the  h  part.  Mostly,  these  two  curves  are  close  to  each  other  and  are  nearly 
indistinguishable.  However,  when  the  algorithm  does  not  converge,  these  curves  can  differ 
widely.  For  analytical  purposes,  there  is  no  essential  difference  in  the  properties  of  these 
two  sequences  of  eigenvalues.  Therefore,  these  curves  need  no  specific  labels. 


Fig.  6.3:  The  illustration  of  the  result  from  Algorithm  4  with  8192  data  points. 


Although  our  program  is  capable  of  generating  the  above  results  within  10  iterations, 
for  any  sample  size  that  is  larger  than  around  213,  our  program  will  require  significantly 
more  time  to  run.  Because  the  use  of  a  smaller  sample  size  merely  implies  a  nonrealistic 
Bdl  which  does  not  compromise  the  characteristic  of  the  CP  that  relates  to  the  channel 
properties,  all  of  our  experiments  were  run  with  sample  size  2U.  We  will  use  the  HT  and 
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TU  channel  profiles  (as  shown  in  Table  6.1)  as  the  multipath  model.  In  order  to  observe  the 
fast  fading  characteristic,  the  channel  needs  to  vary  fast  enough.  At  the  same  time  we  prefer 
to  have  a  positive  sequence  that  closely  approximates  the  original  channel  autocorrelation 
function.  These  together  require  a  large  B 4.  In  fact,  we  use  Bd  =  105, 000.  We  can  obtain 
a  Bd  of  105,000  when  we  travel  at  4,088m/s  while  using  a  carrier  frequency  of  7.7GHz. 
The  parameter  for  frequency  containment  is  set  to  be  /  =  6, 250Hz,  /  =  12, 500Hz,  and 
/  =  25, 000Hz. 
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Chapter  7 

Results  and  Observations 

Some  experimental  results  from  Algorithms  4  and  5  are  presented  in  this  chapter. 
The  multipath  models  as  presented  in  section  6.3  are  used.  The  sequence  of  the  largest 
eigenvalues  obtained  through  iterations  of  the  algorithms  is  used  as  the  measure  of  both 
the  performance  of  the  CPs  and  the  convergence  of  the  algorithms.  We  point  out  that  for 
a  single  path  channel,  the  highest  value  of  the  largest  eigenvalue  from  the  algorithms  is  1. 
This  value  serves  as  a  reference  standard. 

The  results  are  grouped  by  aspects  of  interest  and  are  presented  in  the  following 
subsections:  different  initial  vectors,  frequency  containment,  the  multipath  characteristics, 
the  effect  of  synchronization  uncertainty,  and  the  performance  of  the  signal  sets. 

7.1  Different  Initial  Vectors 

As  we  are  not  certain  about  the  convergence  of  Algorithm  4,  it  is  assuring  to  observe  its 
behavior  for  different  initial  vectors.  The  examples  presented  here  are  with  TU  multipath, 
/3  =  0.3,  /  =  12500Hz,  and  uniformly  distributed  synchronization  delay  with  distribution 
ofW(-4.15  x  10”6,4.15  x  10”6)sec.  As  expected,  Figs.  7.1  and  7.2  show  that  for  different 
initial  vectors,  the  resulting  CPs  are  very  similar  to  each  other  up  to  a  shift  in  time.  The 
largest  eigenvalues  converge  to  0.6994  in  all  cases. 

7.2  Frequency  Containment 

All  plots  in  this  section  use  /  =  12, 500Hz.  They  are  divided  into  four  groups:  Group 
1,  as  shown  in  Figs.  7.3  and  7.4,  comprises  single  path  fading  without  synchronization 
uncertainty;  Group  2,  as  shown  in  Figs.  7.5  and  7.6,  comprises  single  path  fading  with 
synchronization  uncertainty  whose  delay  is  Gaussian-distributed  (a  =  1.124);  Group  3,  as 
shown  in  Figs.  7.7  and  7.8,  comprises  multipath  channel  (HT)  without  synchronization 
uncertainty;  and  Group  4,  as  shown  in  Figs.  7.9  and  7.10,  comprises  multipath  channel 
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(TU)  with  uniformly  distributed  synchronization  delay,  U(— 4.15  x  10  6, 4.15  x  10  6)sec. 
We  can  observe  from  all  these  plots  that  the  signal  and  filter  functions  get  smoother  with 
a  decrease  in  the  largest  eigenvalue  as  8  increases.  We  can  observe  in  Figs.  7.4(b)  and 
7.6(b)  some  numerical  errors  because  the  power  method  does  not  converge  fast  enough. 

7.3  Effect  of  Multipath  Channels 

The  most  interesting  observation  for  multipath  channels  is  that  multiple  paths  may 
increase  the  SNR,  compared  with  single  paths.  All  of  the  plots  in  Figs.  7.11  and  7.12  have 
/3  =  0.  Fig.  7.11(a)  shows  a  plot  with  single  path  fading,  which  has  its  largest  eigenvalue 
converged  to  1.  The  multipath  case  in  Fig.  7.11(b)  has  its  largest  eigenvalue  converged  to 
1.6496. 

We  also  observe  that  the  signals  adapt  to  the  channel  properties.  The  peaks  shown  in 
Fig.  7.11(b)  are  separated  exactly  by  channel  path  delays,  as  expected.  When  there  are 
synchronization  uncertainties,  the  peaky  nature  of  the  pulse  disappears  as  shown  in  Figs. 
7.12(a)  and  (b).  However,  the  largest  eigenvalues  still  have  higher  values  compared  to  that 
of  single  path  settings. 

7.4  Effect  of  Synchronization  Uncertainty 

7.4.1  Fast  Fading 

The  results  in  Figs.  7.13  and  7.14  are  obtained  with  f  =  6250Hz  and  /3  =  0.3  for 
single-path  fading.  The  spread  of  the  signal  pulses  increases  with  the  synchronization 
uncertainty.  At  the  same  time,  the  largest  eigenvalue  or  the  SNR  decreases  as  the  signal 
spread  increases. 

7.4.2  No  Fading 

It  is  interesting  to  see  how  the  synchronization  delay  affects  the  CP.  The  results  can 
be  obtained  by  setting  /3  =  0  and  the  coefficients  of  the  channel  autocorrelation  matrix,  r, 
to  some  constant  values,  which  is  equivalent  to  having  a  constant  autocorrelation  of  the 
channel.  Fig.  7.15(a)  depicts  the  significance  of  the  Match  Filter  Theory.  That  is,  any 


74 


(a)  TU,  Gaussian-distributed  delay  with  a  =  1.124  x  10'  ~6sec 


(b)  TU,  uniform  distribution  with  U(— 4.15  x  10  6,4.15  x  10~6)sec 


Fig.  7.12:  Effect  of  multiple  paths  II. 
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(a)  Perfect  synchronization 


(b)  Gaussian-distributed  delay  with  cr  =  1.124  x  10'  6  sec 


Fig.  7.13:  Effect  of  synchronization  uncertainty  I. 
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(b)  Uniformly  distributed  with  U(— 4.15  x  10  6 , 4.15  x  10  6)sec 

Fig.  7.14:  Effect  of  synchronization  uncertainty  II. 
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signal-filter  pair  will  perform  equally  well  for  a  nonfading-single-path  channel,  as  long  as 
the  signal-filter  pair  is  matched.  On  the  other  hand,  Fig.  7.15(b)  shows  that  if  there  is  a 
synchronization  uncertainty,  the  CP  is  shaped  by  the  distribution  of  delay.  The  effect  of 
synchronization  uncertainty  on  the  shaping  of  CPs  also  occurs  in  multipath  channels  as 
shown  in  Figs.  7.16(b)  and  7.17(b).  Furthermore,  for  the  cases  of  multipath  channels,  a  CP 
is  also  controlled  by  the  path  delays,  even  where  there  is  no  synchronization  uncertainty 
as  shown  in  Figs.  7.16(a)  and  7.17(a). 

7.5  Signal  Sets 

We  discussed  in  section  4.4  that  there  are  two  alternatives  to  apply  the  results  from 
this  dissertation  to  obtain  signal  sets  for  digital  communication.  These  alternatives  are: 
First,  use  the  CPs  generated  with  Algorithm  4  in  a  pulse-position  modulation.  Second, 
use  the  CCPs  generated  with  Algorithm  5.  In  this  section,  we  compare  the  performance 
of  the  signal/filter  pairs  obtained  by  using  our  approach  with  some  traditional  waveforms, 
i.e.,  the  flat-top  pulses  and  the  raised-cosine  pulses. 

7.5.1  CCPs 

Because  of  the  slow  convergence  of  Algorithm  5,  all  plots  in  this  section  are  obtained 
with  512-point-sample  vectors.  We  use  ST  =  1.626  x  10-7s.  This  is  equivalent  to  a  data 
rate  of  12,012  bits  per  second.  In  all  of  the  plots,  we  assume  single  path  fading  channel 
with  perfect  synchronization. 

Figs.  7.18  and  7.19  show  the  change  in  signal/filter  of  CCPs  according  to  different 
values  of  ft.  In  each  of  these  figures,  the  group  of  four  pictures  of  CCPs  is  arranged  into 
two  rows.  On  the  top  row,  the  picture  on  the  left-hand  side  shows  So  while  the  one  on  the 
right-hand  side  shows  ho-  On  the  bottom  row,  the  picture  on  the  left-hand  side  shows  Si 
while  the  one  on  the  right-hand  side  shows  hi .  We  can  observe  from  all  these  plots  that 
the  signal  and  filter  functions  are  getting  smoother  as  /3  increases. 

Fig.  7.20  shows  the  performance  plots  of  probability-of-error  versus  SNR.  On  the 
performance  plots  of  joint  CPs,  we  also  show  that  of  binary  flat-top  pulse  and  raised- 
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(a)  Single  path,  perfect  synchronization 


(b)  Single  path,  uniformly  distributed  delay  with  U(— 4.15  x  10  6 , 4.15  x 
10_6)sec 


Fig.  7.15:  CP  for  channels  without  fading  I. 


300 


(a)  TU,  perfect  synchronization 


« 
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(b)  TU,  uniformly  distributed  delay  with  U(— 4.15  x  10  6 , 4.15  x  10  6)sec 


Fig.  7.16:  CP  for  channels  without  fading  II. 
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SNR  vs  P(e) 


Fig.  7.20:  SNR  vs.  P(e),  CCPs  for  single  path,  perfect  synchronization. 

cosine  pulse  in  pulse  position  modulation  (PPM)  pairs  with  a  support  of  8.35  x  10-5s  in 
the  same  channel  for  comparison.  The  observations  are:  First,  for  CCPs,  the  smaller  the 
/?,  the  better  the  performance.  Second,  there  are  error  floors  for  the  CCPs.  Third,  before 
reaching  their  error  floor,  CCPs  with  f3  =  0.1  and  /3  =  0.3  have  a  performance  gain  of  2  dB 
to  3  dB  over  the  raised-cosine  pulse,  and  have  a  8  dB  to  9  dB  gain  over  the  flat-top  pulse. 
CCPs  with  (3  =  0.5  have  a  performance  marginally  better  than  that  of  the  raised-cosine 
pulse,  while  that  with  /3  —  0.7  have  a  performance  loss  of  around  2  dB  over  that  of  the 
raised-cosine  pulse. 

From  all  of  the  experiments  we  did  for  different  fading  channels,  there  is  a  similar  trend 
of  relationships  between  the  performance  of  CCPs,  raised-cosine  pulse  and  flat-top  pulse 
as  observed  above.  Another  observation  is  different  initial  vectors  can  generate  CCPs  with 
very  different  performance.  We  used  four  different  initial  vectors  for  each  set  of  channel 
parameters  and  ran  Algorithm  5  for  60  iterations  with  each  of  these  initial  vectors.  The 
best  CCPs  obtained  with  each  of  the  different  initial  vector  can  differ  by  more  than  ten 
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f„,ds  in  terms  of  probability-error  for  the  same  set  of  channei  parameters.  Considering 
the  small  number  of  experiments  we  did,  the  possible  best  CCPs  may  perform  far  better 

than  these  CCPs  we  found. 

7.5.2  PPM  with  CPs 

Results  presented  in  Figs.  7.21  through  7.24  are  obtained  with  f  =  12,500Hs  an 

ST  =  1 626  x  10-7s.  From  all  the  plots,  we  observe:  Erst,  a  5  to  10  dB  gam  in  performance 
for  the  raised-cosine  pulse  over  the  ttat-top  pulse  for  different  Ming  and  synchronisatron 
situations;  and  second,  CPs  with  V  =  0.1  through  „  =  0.5  have  very  similar  performance 
while  the  one  with  Mu  -  0.7  is  generally  inferior  to  the  other  CPs.  Figure  7.21  shows 

the  performance  of  the  CPs  for  TU  multipath  fading  model  wrth  a  umformly  distn  u  e 

,  .  jj(_ 4  15  x  io~6, 4.15  x  10-6)sec.  The  waveforms  of  these 

synchronization  delay  such  that  ( 

.  .  ,  7  aIld  7  8.  There  is  an  8  to  9  dB  gain  in  performance  of  all  these 

CPs  are  shown  in  ! igs.  <•<  ana  1  -°- 

CPs  over  the  raised-cosine  pulse.  i  .  . 

,  r  nt  thp  CPs  for  HT  multipath  fading  model  without 

Figure  7.22  shows  the  performance  of  the  <^s  tor  n  v 

i  •  7  7  and  7  8  We  also  observe 

synchronization  delay  whose  waveforms  are  shown  in  g  •  •  ' 

•  •  r  anrp  nf  the  CPs  over  the  raised-cosine  pulse.  Now,  t  ere  is 

a  10  to  12  dB  gam  m  performance  of  the  o 

less  than  2  dB  difference  between  the  CP  with  (3  =  0.7  and  the  other  CPs. 

Figure  7.23  shows  the  performance  of  the  CPs  for  single  path  fading  with  a  Gauss, an- 
distributed  synchronisation  delay  (,  =  UM),  whose  waveforms  for  different  f.  are  shown 
ln  Figs  7.5  and  7.6.  In  this  figure,  we  observe  that  the  flat-top  pulse  reaches  its  error  floor 
at  the  SNR  of  65  dB.  We  also  observe  a  10  to  13  dB  gain  in  performance  of  the  CPs  over 
the  raised-cosine  pulse.  This  time,  there  is  a  3  dB  difference  between  the  CP  with  f  =  • 

and  the  other  CPs. 

Figure  7.24  shows  the  performance  of  the  CPs  for  single  path  fading  with  a  per  ec 
synchronization  whose  waveforms  for  different  values  of  f  are  shown  in  Figs.  7.3  and  7.4. 
we  observe  a  10  to  15  dB  gain  in  performance  of  the  CPs  over  the  raised-cosine  pulse, 
while  the  CP  with  y  0.7  has  a  5  dB  difference  from  the  other  CPs. 
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7  6  The  Convergence  of  Algorithms 

'  In  s„me  cases,  Algorithms  4  and  5  do  not  converge.  Among  all  unknown  causes  of  the 
failure  for  these  algorithms  to  converge,  there  are  two  major  possible  causes.  Fustly,  the 
power  method  or  the  conjugate  gradient  method  may  not  have  converged  raprd  enough  to 
the  correct  eigenvalue/eigenvector  by  the  time  the  preset  number  of  iterations  has  been 
exhausted.  This  situation  may  occur  if  the  largest  eigenvalues  are  too  close  together  for  the 
power  method,  or  the  smallest  eigenvalues  are  too  close  together  for  the  conjugate  gradrent 
method  (see  discussion  in  section  6.4).  Secondly,  the  averaging  effect  of  synchromsat.on 


uncertainty  counteracts  with  the  convergence. 
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Chapter  8 

Conclusion  and  Discussion 

m  this  chapter,  the  contributions  made  in  this  dissertation  are  summarized.  Directions 
for  future  study  are  discussed. 

8.1  Summary 

In  this  communication  theory  research,  we  were  able  to  characterize  the  optimal  detec¬ 
tor  and  signal  by  maximizing  the  SNR  tor  communication  through  multiplicative  channels, 
which  may  be  slow  or  fast  fading.  In  addition,  this  characterisation  is  general  enough  to 
include  single-path  and  multipath  fading  with  different  synchronisation-delay  profiles. 

The  design  criterion  of  maximizing  SNR  with  feed  energy  leads  to  a  constram  op¬ 
timisation  problem.  The  optimization  was  performed  by  using  the  Lagrange  mult.phers 
in  a  function  space.  We  were  able  to  characterize  the  optimal  detector  as  the  ergen- 
function  corresponding  to  the  largest  eigenvalue  of  an  operator  with  kernel  bemg  the 
product  of  the  channel  autocorrelation  function  and  a  transformation  of  the  transmu¬ 
ted  signal.  The  transformation  of  the  transmitted  signal  depends  on  the  multipath  and 

the  synchronization-delay  profiles. 

We  then  derived  the  optimal  signal  and  detector  as  a  pair,  this  time  optimising  the 
variables  simultaneously  using  the  Lagrange  multipliers.  The  conclusion  we  obtamed 
was  that  the  signal  and  the  filter  must  be  similar  to  each  other  in  amplitude,  but  they 
ought  differ  in  sign.  However,  the  resulting  expression  did  not  provide  a  way  to  determine 
the  optimal  functions  physically  besides  a  list  of  properties  as  in  section  3.4. 

Based  on  the  symmetric  role  of  the  filter  and  signal  in  a  communication  system, 
we  derived  an  iterative  algorithm  (Algorithm  1)  which  allowed  -  to  generate  an  opt.mal 
consonant  pair  from  random  seeds.  This  iterative  algorithm  contains  two  alternating  steps. 
!  an  optimal  detector  to  generate  an  updated  filter  function  for  the  previously  generated 
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for  the  previously  generated  filter  function. 

For  single-path  fading,  the  above  iterative  algorithm  always  generates  consonant  pairs 

that  approach  delta  distributions.  Physically  this  is  very  reasonable.  If  there  is  an  uncer¬ 
tainty  in  the  channel  status,  a  transmitted  signal  should  use  as  little  time  of  the  dranne 
as  possible.  However,  a  taction  that  is  d.  to  the  delta  distribution  is  not  pract.cal,  be¬ 
cause  it  requires  infinite  band  width  to  transmit.  So,  we  introduce  a  frequency  contamment 

constraint  to  Algorithm  1  and  derived  Algorithm  2. 

In  order  to  communicate  over  a  fast  fading  channel,  the  use  of  amplitude  and  phase 
modulation  is  inappropriate.  Something  akin  to  orthogonal  signaling  is  necessary.  A 
procedure,  Algorithm  3  was  developed  to  generate  two  CCPs  such  that  the  correspond, ug 
functions  were  neariy  orthogonal  among  the  couiugate  pairs.  Algorithm  3  required  two 
copies  of  Aigorithm  2  being  executed  in  ptallel  with  the  additional  requirement  that  the 
updated  consonant  pair  from  one  execution  was  related  to  the  eigenfunction  corresponding 

to  the  smallest  eigenvalue  of  the  other  execution. 

Furthermore,  we  derived  Algorithms  4  and  5,  the  discrete  version  of  Algorithms  2  and 

3  on  which  we  did  all  the  experiments. 

’  to  help  us  compare  CCPs  generated  by  our  procedure  with  other  standard  signaling 
waveforms,  we  aUo  derived  the  density  and  distribution  tactions  for  quadratic  forms 
of  complex  normal  vectors.  The  equation  we  obtained  can  be  evaluated  by  usmg  the 
eigenvalues  of  a  single  matrix.  This  matrix  is  the  product  of  the  correlation  matnx  of  the 
normal  vector  and  the  Hermitian  of  the  quadratic  form.  This  is  an  extens.on  of  prevous 
results  as  most  of  these  resuits  would  require  us  to  evaluate  some  infinite  sequence  or 
only  applicable  to  some  specific  values.  A  partition  of  unity  found  in  the  evaluat.ou  of 
the  distribution  function  at  point  zero  is  useful  in  simplifying  some  expressions  used  m 
standard  digital  communication  reference  books. 
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8.2  Future  Work 

This  dissertation  provides  a  new  approach  to  signal  and  detector  design  for  commu¬ 
nication  through  fading  channels.  There  are  many  directions  that  can  be  taken  from  this 
point  both  theoretically  and  practically.  Here  we  list  a  few  of  them. 

1.  Convergence  of  the  algorithm:  It  is  easy  to  show  the  algorithm  for  single-path  fading 
channels  when  there  is  no  synchronization  uncertainty  with  /3  =  0  or  /3  =  1  has  a 
monotonic  increasing  largest  eigenvalue.  The  convergence  of  the  algorithm  for  the 
above  and  other  situations  remains  unknown.  Knowing  the  limit  of  the  convergence 
of  the  algorithm  will  help  to  extend  the  analysis  by  applying  similar  technique  to 
more  the  complex  models. 

2.  Quadratic  forms:  The  development  of  the  expression  of  the  density  and  the  distribu¬ 
tion  functions  of  quadratic  forms  with  zero  mean  normal  vectors  is  sufficient  for  this 
dissertation.  However,  to  make  the  expression  for  the  density  and  the  distribution 
functions  truly  useful,  e.g.,  applicable  to  the  slow  fading  situations,  an  extension  of 
these  equations  to  include  nonzero  mean  normal  vectors  is  necessary  and  is  highly 
valued. 

3.  Error  floor:  The  phenomenon  of  error  floor  is  discussed  in  the  literature,  including 
those  written  by  Hansson  [108]  and  Simon  [92].  The  concept  of  error  floor  is  referred 
to  as  “irreducible  bit  error  probability”  or  “error  floor.”  The  error  floor  is  the  limiting 
behavior  of  P(e)  >  0  as  the  SNR  approaches  infinity  for  some  signaling  schemes  with 
their  corresponding  receivers.  Simon  pointed  out  that  for  M-ary  Differential  phase- 
shift-keying,  the  error  floor  exists  in  any  channel  with  fading  autocorrelation  that 
is  not  equal  to  one  (fast  fading).  Hansson,  on  the  other  hand,  showed  that  for 
orthogonal  signal/filter  pairs,  the  error  floor  does  not  exists.  Results  from  section 
5.5  can  be  used  to  extend  the  characterization  of  the  signal/filter  pairs  that  do  not 


have  an  error  floor. 
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4  Consonant  pairs  for  slow/nonfading  channels:  We  observed  from  the  experinrents 
that  there  are  relationships  between  the  synchronisation  uncertainty  and  the  shapes 
of  a  consonant  pair.  This  is  an  area  that  may  have  significant  practical  unportance, 
as  synchronization  uncertainty  is  unavoidable  in  communication. 


5  Analysis  of  the  consonant  pair:  We  can  observe  that  consonant  pans  of  ddferent 
Channels  share  many  similar  properties.  If  we  can  analyze  and  successfully  represent 
the  consonant  pairs  of  most  common  channels  with  some  specihc  family  of  functrons, 
e  g  wavelets,  we  have  a  flexible  communication  system  that  can  easily  adapt  to  these 


common  channels. 
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Appendix  A 

Derivative  of  a  Functional  with  Respect  to  a  Function 

Let  function  u  G  £2  such  that  u  :  E  -»  R,  and  let  functional  F  :  £2  ->  R  such  that 

Ffu)  G  R  Then  G  £2,  the  derivative  of  jP  with  respect  to  u  at  uq  G  £2  can  be 

v  '  du 

defined  as  a  linearization  of  F  at  «o  in  Gateaux  sense  [123]  as  follows: 

F(u0  +  eU)  =F(u  0)  +  eU)  +  o(e) 

=f(«o)  +  e@,t0  +  o(4 

du 

where  e  G  M,  and  U  G  £2  is  any  function,  and 

fiEp. £,u)=  [  (^r^-)W(t)  dt. 

du  J  du 

As  e  -4  0,  o(e)  ->  0,  from  (A.l)  and  (A.2),  we  have 

[ dt  =  lim  f  ("°  ~  (V)  ~  F{n) .  (A.3) 

J  v  du  w  w  e->o  e 

Assume  F(u)  is  differentiable  at  u0  and  F{u0)  is  a  local  maximum  of  F.  Then  for 
e  >  0  and  when  e  is  small,  there  exists  K  G  M  such  that 

F(„„  -  eCQ  -  FM  <R  f  ,muo)){t)m  dt  <  F(»q  -  rip  -  FM  (A  4) 
— e  J  du  € 

As  (A.4)  is  true  for  all  U  €  £2,  especially  for  U  =  0,  we  have 

/*  dF(uo 2)(t)l/(t)  =  0.  (A.5) 

J  du 


(A.l) 


(A.2) 


Additionally,  (A.5)  implies 

^(«0)  _ n 

du 

since  it  is  true  for  all  U  €  £2-  The  same  argument  with  e  <  0  can  show  that 
for  F(u0)  being  a  local  minimum. 


(A.6) 


dF(uo) 


=  0 


du 
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The  definition  of  derivative  of  a  functional  with  respect  to  a  function  as  in  (A. 3)  can 
be  extended  to  complex  functions.  Let  h  G  £2  such  that  h  :  R  — t  C,  and  functional 
p  .  £2  £  such  that  F(h)  G  R.  The  extension  may  be  done  in  the  following  manner.  Let 

h(t)  =  u(t)  +  jv(t),  where  both  u(t)  and  v(t)  are  real,  we  have  1 


dFjh)  _  1  (dFjk)  _  .9F{h)\ 


dh 


2  \  du 


dv  )  ' 


and 


where 


dF(h)  _  1  (dF{h)  ,  .dF(h)\ 


dh*  2  V  d 


u 


+  3- 


dv  )  5 


(A.7) 


(A. 8) 


dF(hl  and  dF^  are  computed  as,  for  e  G  M,  and  U,V  £  £2  being  real  valued, 


d 


u 


dv 


ho 


,v) = dt 


and 


ho 


F{h0  +  eU)  —  F{h$) 

=  lim - - - 

c — »0  € 

F({uq  +  eU)  +  jv 0)  -  F(ho) 

=  lim  - - j 

c— ^0  e 


.V)  =/(^2i)(()F(t)  <i( 

F(h0+ jeV)  -  F(ho) 

=  lim - - - 

€ — >0  6 

,.  ^(uo  +j{vo  +  eV))  -F(ho) 

—  hm  — - — • 

e-fO  e 


(A.9) 


(A. 10) 


To 


demonstrate  the  above  computation  with  an  example,  suppose  H  ( h )  —  /  h{t)h*{t)  dt , 


then 


if  ((tio  +  €l/)  +  3V 0)  ~  H(ho) 

=  j [(uo{t)  +  eU{t))  +  J«o(t)][(uo(t)  +  eCf(t))  -  j«o(i)]  dt  -  H{h0) 
=  J {{u0{t)  +  eU(t))2  +  vg(t)>  dt-  J {u20{t)  +  v20{t)}  dt 

«2e  J  uo(t)U(t)  dt, 


(A.11) 


1  Following  Haykin  [120],  whom  defines  in  this  manner  for  the  complex  derivatives  with  discrete  vectors, 
which  in  turn  follows  Schwartz  [124]. 
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so, 


dFjhp ) 
du 


)(t)  =  2u0(t). 


(A.12) 


Also, 


H(u0  +  j(vo  +  eV))  —  H(ho) 

=  J [uoW  +  jM<)  +  eV(i))][it0(i)  ~  j(v o{t)  +  eV{t))}  dt  -  H{h0) 
=  f  {«oW  +  MO  +  *V{t))2}  dt-  j {u20(t)  +  wg(t)}  dt 

w2e  J  vQ(t)V(t)  dt , 


(A.13) 


gives 


dFjhp) 

dv 


)(t)  -  2 v0{t). 


(A.14) 


dF{h) 

dh 


(t)  =  h*0(t) 

ho 


and 


dF{h) 

dh* 


(t)  =  ho{t). 

ho 


(A.15) 


That  is, 
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Appendix  B 

Inverse  of  a  Symmetric  Kernel 

Let  k( X,r)  =  l[g{X,T)},  where  !(■)  stands  for  the  inverse  of  the  operator  g(X,r).  So 
that  we  have 

J  k(X,T)g{T,u)d,T  =  8(X  -  u).  (B.l) 

For  the  case  g{ X,  r)  =  g{ X  -  r )  and  k{ X,  r)  =  k{ X  -  r),  (B.l)  becomes 

J  k(t  —  u)g(u)du  =  S(t)  (B-2) 

where  u  =  X  -  r.  Observe  that  the  left-hand  side  of  (B.2)  represents  the  convolution 
between  k(t)  and  g(t).  By  performing  the  Fourier  transform  to  both  left  hand  side  and 
right  hand  side  of  (B.2),  we  obtain 

K(u)G{u)  =  1,  (B.3) 

or 

<R4) 

where  K(u)  and  G{oj)  are  the  Fourier  transforms  of  k(t)  and  g(t),  respectively.  If  G{ui)  =  0 
for  some  uj  and  is  bounded  below,  we  use  K(iS)  =  ^  the  approximation  of  the  inverse 

operator  kernel,  with  I  being  some  constants  such  that  G{oj)  +  I  >  0. 


Appendix  C 
Programs 


Matlab  Codes 


Algorithm  4 

function  [si,SH,EV,f lag]  =  alg(signal_length,deltaT,Bd, beta, xi, omega, K, mu, 
'/,  [si,SH,EV,flag]  =  alg(signal_length,deltaT,Bd, beta, xi, omega, K, mu, N) 

7. 

•/,  This  program  implements  Algorithm  2  in  the  dissertation  which  having 
7,  Algorithm  1  as  a  special  case  when  beta  =  0 

7. 

7,  Output: 

7. 

7,  Si:  initial  signal 
7.  SH:  final  S,  and  H  pair 

7.  EV:  eigenvalue  recorded  through  out  the  process 
7,  flag:  flag  =  1  when  there  is  any  decrease  in  ev 

7. 

7. 

7.  Input: 

7. 

7,  signal_length:  number  of  points  in  the  signal. 

7,  deltaT  :  sampling  time 

7,  Bd:  Doppler  freq 

7,  beta:  weight  number  [0,1] 

7,  xi  =  [  phi_l  xi_l] 

7,  [  phi_2  xi_2] 

7.  [  phi_3  xi_3] 

7,  [  phi_4  xi_4] 

7,  : 

•/,  phi_i  is  the  delay,  and  xi_i  is  the  channel  gain  for  path  i 

7,  omega:  bandwidth  in  angular  speed 

7,  K:  span  of  the  random  variable 

7,  mu:  synchronization  delay 

7.  0 ,  delta 

7,  1,  uniform 

7,  2,  Gaussian 

7,  N:  number  of  iterations 

7. 

sxi  =  xi; 

smxi  =  max(sxi ( : , 1) ) ; 
hxi=[-l*xi(: ,1)  xi(:,2)]; 
hmxi  =  min(hxi(:  ,D) ; 
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EV  =  □  ; 
oldev  =0; 
flag  =  0; 

SH  =[] ; 

r  =  Cal_r  (bet a , Bd , s ignal_length , delt aT)  ; '/,  autocorrelation  function  of  the  channel 
rl  =  Cal_r(0,Bd,signal_length, deltaT) ; 

W  =  Cal_W  (beta ,  omega,  signal  .length,  deltaT) ;  '/.  2beta  w  sinc(2pi  w(tau-lambda)) 

WO  =  zeros(size(W)) ; 
h  =  zeros (signal_length, 1) ; 

Z1  =  floor (signal_length/4) ; 

Z  =  zeros (Zl, 1) ; 

s  =  [Z  ;randn(  (signal_length“2*Zl)  ,  1) ;  Z]  ;  */,  leave  room  for  SPmu  and  shifts 
s  =  normal_v(s , deltaT)  ; 
si  =  s; 
for  n  =  1:N, 

fprintf  ( ’XnNumber  of  iteration :  5£d\n*  ,n)  ; 
s  =  SPmu(s,K,mu) ; 

[h,hev]  =  Maxeig_p _method (r , s , sxi , smxi , W , 0) ; 
h  =  normal_v(h, deltaT) ; 
hev  =  h*  *Rxb(rl , s ,sxi ,smxi ,W0,0 ,h) ; 
if  hev  <  oldev, 
flag  =  1 
end 

oldev  =  hev; 
h  =  SPmn(h,K,mu) ; 

[s ,sev]  =  Maxeig_p.method(r,h,hxi,hmxi,W,0) ; 
s  =  normal^v (s, deltaT) ; 
sev  =  s,*Rxb(rl,h,hxi,hmxi,W0,0,s) ; 
if  sev  <  oldev, 
flag  =  1 
end 

oldev  =  sev; 

EV  =  [EV;  sev  hev]; 

SH  =  [s  h]  ; 

end 

EV  =  EV*deltaT*deltaT ; 

Algorithm  4  for  Non-Fading  Channel 
function  [si,SH,EV,f lag]  = 

alg.slow (signal.length , deltaT , Bd , beta , xi , omega ,K,mu,N) 

V,  [si ,SH,EV,f lag]  -  alg_slow (signal.length , deltaT ,Bd, beta, xi , omega, K , mu ,N) 

•/. 

*/,  This  program  implements  Algorithm  2  for  the  case  when  there  is  no  fading. 

y. 

'/,  Output: 

% 

*/,  Si:  initial  signal 
*/,  SH:  final  S,  and  H  pair 

'/,  EV :  eigenvalue  recorded  through  out  the  process 
*/,  flag:  flag  =  1  when  there  is  any  decrease  in  ev 


108 


7. 

% 

7*  Input: 

7. 

*/,  signal.length:  number  of  points  in  the  signal. 

7,  deltaT  :  sampling  time 

7.  Bd:  Doppler  freq 

7,  beta:  weight  number  [0,1] 

7,  xi  =  [  phi.l  xi_l] 

7#  [  phi_2  xi_2] 

7,  [  phi. 3  xi_3] 

7,  [  phi_4  xi_4] 

7.  : 

•/,  phi.i  is  the  delay,  and  xi.i  is  the  channel  gain  for  path  i 

7.  omega:  bandwidth  in  angular  speed 

7,  K:  span  of  the  random  variable 

7,  mu:  synchronization  delay 

7.  0,  delta 

7,  1,  uniform 

7,  2,  Gaussian 

7*  N:  number  of  iterations 

7. 

sxi  =  xi; 

smxi  =  max(sxi( : , 1) ) ; 
hxi=[-l*xi(: ,1)  xi(:,2)]; 
hmxi  =  min(hxi(: ,1)) ; 

EV  =  []  ; 
oldev  =0; 
flag  =  0; 

SH  =[]; 

r  =  ones(signal_length,l) ; 

W0  =  zeros (signal.length, 1) ; 
h  =  zeros (signal.length, 1) ; 

Z1  =  floor (signal_length/4) ; 

Z  =  zeros(Zl,l) ; 

s  =  [Z  ;randn(  (signal_length-2*Zl)  ,  1)  ;Z]  ;  7,  leave  room  for  SPmu  and  shifts 
s  =  normal_v (s, deltaT) ; 
si  =  s; 
for  n  =  1:N, 

f printf  ( * \nNumber  of  iteration : 7.d\n>  ,n)  ; 
s  =  SPmu(s,K,mu) ; 

[h,hev]  =  Maxeig_p_method(r , s , sxi , smxi , W0, 0) ; 
h  =  normal. v(h, deltaT) ; 
hev  =  h } *Rxb (r , s , sxi , smxi , W0 , 0 , h) ; 
if  hev  <  oldev, 
flag  =  1 
end 

oldev  =  hev; 
h  =  SPmu(h,K,mu) ; 

[s,sev]  =  Maxeig_p_method(r ,h,hxi ,hmxi ,W0,0) ; 
s  =  normal. v(s , deltaT) ; 
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sev  =  s’*Rxb(r,h,lixi,hinxi,WO,0,s) ; 
if  sev  <  oldev, 
flag  =  1 
end 

oldev  =  sev; 

EV  =  [EV;  sev  hev]  ; 

SH  =  [s  h]  ; 

end 

EV  =  EV*deltaT*deltaT ; 


Algorithm  5 
function  [MaxSH,pe30]  = 

algpairpe30(signal_length, deltaT,Bd, beta, xi, omega, K, mu, N1,N2) 

*/,  [MaxSH,pe30]  =  algpairpe30(signal_length, deltaT, Bd, beta, xi, omega, K,mu,Nl,N2) 

•/. 

*/. 

'/,  Output: 

•/. 

V. 

*/. 

'/,  Input: 

•/. 

*/.  signal_lengtb:  number  of  points  in  the  signal. 

'/,  deltaT  :  sampling  time 
*/.  beta:  weight  number  [0,1] 

'/,  xi  =  [  phi_l  xi_l] 

*/,  [  phi_2  xi_2] 

*/,  [  phi_3  xi_3] 

*/,  [  phi_4  xi_4] 

*/.  : 

yt  phi_i  is  the  delay,  and  xi_i  is  the  channel  gain  for  path  i 

*/,  omega:  bandwidth  in  angular  speed 

'/,  K:  span  of  the  random  variable 

'/,  mu:  synchronization  delay 

*/,  0,  delta 

'/,  1,  uniform 

*/,  2,  Gaussian 

'/,  Nl:  number  of  time  using  random  initial  vector 
'/,  N2:  number  of  iterations  for  each  initial  vector 

sxi  =  xi; 

smxi  =  max(sxi(: ,1)) ; 
hxi=[-l*xi(: , 1)  xi(:,2)]; 
hmxi  =  min(hxi(:  ,1))  ; 

EV  =  □  ; 
oldev  =0; 
flag  =  0; 

SH  =[] ; 

r  =  Cal_r (beta, Bd,signal_length, deltaT)  ; '/.autocorrelation  function  of  the  channel 
rO  =  Cal_r(0, Bd,signal_length, deltaT) ; 


no 


W  =  Cal_W (beta , omega , signal_length , deltaT) ;  */, 2bw  sinc(2w(tau-lambda)) 

WO  =  zeros (size (W) ) ; 
if  beta  ==  0 

IW  =  W;  ’/.doesn't  matter 

else 

IW  =  -9 .091*W/beta; 

IW(1)  =  IW(1)  +  61500620;  */.  10* (1 . 626e-7) ~ (-1) ; 

IW  -  beta*IW; 
end 

for  n  =  1:N1, 

f printf ( ' \nNumber  of  iteration t’/.dXn' ,n) ; 
s  =  zeros (signal_length, 1) ; 

Z  =  zeros (Zl,l) ; 

h  =  randn((signal_length) ,  1)  ;  ’/. 
h  =  normally (h, deltaT) ; 
h  =  SPmu(h,K,mn) ; 

[s0,e01]  =  Maxeig_p_method(r ,h ,hxi ,hmxi ,W,0) ; 
sO  =  normal_v(s0, deltaT) ; 

[si ,e02]  =  Mineig_cg_method(r ,h ,hxi ,hmxi , IW , 10) ; 
si  =  normally (s 1 , deltaT) ; 
mperr  =  1; 

for  n2  -  1:N2 

f pr intf ( > \b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b ’ ) ; 
f printf  ( ’ \b\b\b\b\b\b\b\b\b\b\b\b\b\b\bNl : */,d  N2:*/,d  Hpe:*/,f ’ ,n,n2, mperr) ; 
s  =  SPmu(sO,K,mu) ; 

[hOO , ebOmax]  =  Maxeig_p_method (r ,  s , sxi , smxi , W , 0) ; 

[hl0,eh0min]  =  Mineig_cg_method(r ,s , sxi, smxi, IW, 10) ; 

s  =  SPmu(sl,K,mu) ; 

[hll ,ehlmax]  =  Maxeig_p_method(r ,s, sxi, smxi ,W,0) ; 

[h01,ehlmin]  =  Mineig_cg_method(r,s, sxi, smxi, IW, 10) ; 

h0_new  =  hOO  +  hOl; 
hl_new  =  hlO  +  hll; 
hO  =  normal_v (hO_new, deltaT) ; 
hi  =  normal.v (hl_new, deltaT) ; 

h  =  SPmu(hO ,K,mu) ; 

[sOO ,esOmax]  =  Maxeig_p_method(r ,h,hxi,hmxi,W,0) ; 

[slO,esOmin]  =  Mineig_cg_method(r ,h,hxi,hmxi, IW, 10) ; 

h  =  SPmu(hl,K,mu) ; 

[sll ,eslmax]  =  Maxeig_p_method(r ,h,hxi,hmxi,W,0) ; 

[s01,eslmin]  =  Mineig„cg_method(r ,h,hxi,hmxi,IW, 10) ; 

sO_new  =  sOO  +  sOl; 
sl_new  =  slO  +  sll; 
sO  =  normally (sO_new, deltaT) ; 
si  =  normal_v(sl_new, deltaT) ; 
sh  =  [sO  si  hO  hi] ; 


'/,  keep  the  one  with  smaller  P(e)  at  30dB 
nperr  =  perr(sh,r0,sxi,smxi,W0,K,mu,30,deltaT) ; 
if  mperr  >  nperr 
MaxSH( : , : ,n)  =  sh; 
pe30(n)  =  nperr; 
mperr  =  nperr; 
end 
end 

end 

P(e)  for  Quadratic  Receiver 

function  [err]  =  perr(SH, r, xiphi, mxiphi, W,K, mu, SNR, deltaT) 

'/,  [err]  =  perrpair(SH,R, XiPhi, SNR) 

*/.  SH  =  [sO  si  hO  hi] 

'/,  xiphi  =  [  phi_l  xi_l] 

'/,  [  phi_2  xi_2] 

V,  [  phi_3  xi_3] 

7.  [  phi_4  xi_4] 

*/.  : 

y§  phi_i  is  the  delay,  and  xi_i  is  the  channel  gain  for  path  i 

’/,  mxiphi:  max  delay 

DN  =  2; 

CO  =  zeros(2,2) ; 

Cl  =  zeros(2,2) ; 

MO  =  zeros (2, 2) ; 

Ml  =  zeros(2,2) ; 

H  =  zeros (2, 2); 

Q  =  zeros(2,2) ; 
sO  =  normal_v(SH(: ,1) ,deltaT) ; 
si  =  normal_v(SH( : ,2) ,deltaT) ; 
hO  =  normal_v(SH( : ,3) ,deltaT) ; 
hi  =  normal_v(SH(: ,4) ,deltaT) ; 
sO  =  sO*sign(sO’*hO) ; 
si  =  sl*sign(sl’*hl) ; 
s  =  SPmu(sO,K,mu) ; 

MO  =  [hO ’*Rxb(r,s, xiphi, mxiphi ,W,0,h0)  hO’*Rxb(r, s, xiphi, mxiphi, W,0, hi) ; 

hi ’*Rxb(r,s, xiphi, mxiphi, W,0,h0)  hi ’ *Rxb(r,s, xiphi, mxiphi, W,0, hi)] 
MO  =  M0*deltaT"2; 
s  =  SPmu(sl,K,mu) ; 

Ml  =  [hO’*Rxb(r,s, xiphi, mxiphi, W,0,h0)  hO’*Rxb(r,s, xiphi, mxiphi, W,0, hi) ; 

hi’ *Rxb(r,s, xiphi, mxiphi, W,0,h0)  hi’ *Rxb(r,s, xiphi, mxiphi, W,0, hi)] 
Ml  =  Ml*deltaT“2; 

H  =  [1  hO’*hl*deltaT;  hl’*hO*deltaT  1]; 

PO  =  0; 

PI  =  0; 

N2  =  10“  (-SNR/ 10)  ;  ‘/.N/2 
CO  =  MO  +  N2*H; 

Cl  =  Ml  +  N2*H; 

xi=  log(det(C0)/det(Cl)); 


Q  =  inv(Cl)-inv(CO) ; 
a  =  eig(CO*Q) ; 
b  =  eig(Cl*Q) ; 
if  xi  <=  0, 
for  I  -  1:DN, 
temp  =  1 ; 
if  a(I)  <  0, 
for  n  =  1:DN, 
if  n  “=  I, 

temp  —  temp*(l/ (l-(a(n)/a(I)))) , 
end 
end 

temp  =  temp*exp(-xi/a(I)) ; 

P0  =  P0  +  temp; 
end 

temp  =  1; 
if  b(I)  <  0, 
for  n  =  1:DN, 
if  n  ~=  I, 

temp  =  temp*(l/ (l-(b(n) /b(I)))) ; 
end 
end 

temp  =  temp*exp(-xi/b(I)) ; 

PI  =  PI  +  temp; 
end 
end 

err  =  0.5*(1+P0-P1) ; 
else 

for  I  =  1 :DN, 
temp  =  1; 
if  a(I)  >  0, 
for  n  =  1:DN, 
if  n  "=  I, 

temp  =  temp*(l/(l-(a(n)/a(I)))) ; 
end 
end 

temp  =  temp*expC-xi/a(I)) ; 

P0  =  P0  +  temp; 
end 

temp  -  1 ; 
if  b(I)  >  0, 
for  n  =  1:DN, 
if  n  "=  I, 

temp  =  temp*(l/(l~(b(n)/b(I)))) ; 
end 
end 

temp  =  temp*exp(-xi/b(I)) ; 

PI  =  PI  +  temp; 
end 
end 

err  =  0.5*(1-P0+P1) ; 
end 


Power  Method 


function  [V,ev]  =  Maxeig_p_method(r_alpha,s,xi,mxi,W,kk) 

*/.  [V ,  ev]  =  Maxeig.p .method (r.alpha, s ,xi  ,mxi ,W,kk) 

7.  .  .  _ 

'/,  This  function  returns  the  largest  eigenvalue  with  the  corresponding 
*/,  eigenvector  using  power  method 

7. 

7,  r.alpha:  R.Xalpha 
'/,  xi  =  [  phi.l  xi.l] 

•/,  [  phi_2  xi_2] 

*/,  [  phi. 3  xi_3] 

*/,  [  phi_4  xi_4] 

•/'  phi.i  is  the  delay,  and  xi.i  is  the  channel  gain  for  path  i 

'/,  mxi:  the  largest  delay 

*/.  W  :  2\beta\omega\sinc(2\pi\omega\delta  T  n) 

*/,  kk  :  a  constant  that  allow  us  to  manipulate  the  eigenvalues  of  the 
'/,  kernel 

7. 


V  =  randn(size (r.alpha) ) ; 

V  =  V/norm(V) ; 

V.old  =  V; 
arg  =  1; 

i  =  0; 

while  (arg  >  le-6)  &  (i<50) 
i  =  i+1; 

V  =  Rxb (r .alpha, s,xi, mxi, W, kk, V) ; 

V  =  V/norm(V) ; 

arg  =  norm(V-V.old) ; 

f pr intf ( ’ \b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b ’ ) ; 
f printf ( * \b\b\b\b\b\b\b\bMaxeig.p_method : ^d  ^f  ,i,arg)j 
V.old  =  V; 
end 

b  =  Rxb(r_alpha,s,xi,mxi,W,kk,V) ; 
ev  =  norm(b)  -  kk; 


Conjugate  Gradient  Method 

function  [V,ev]  =  Mineig_cg_method(r_alpha,s,xi,mxi,W,kk) 

7.  [V , ev]  =  Mine ig.cg.method (r.alpha ,s,xi,mxi,W,kk) 

7,  Using  Conjugate  Gradient  Method.  This  method  applies  to  symmetric 
7.  matrix  R  which  solves  Rx  =  b. 

7,  P.311  of  James  W.  Demmel  "Applied  numerical  linear  algebra" 

7. 

7,  This  program  returns  the  smallest  eigenvalue  and  the  corresponding 
7,  eigenvector  of  R 

7. 

7,  r.alpha:  R.Xalpha 
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7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 

7. 


xi  =  [  phi-1  xi_l] 
[  phi_2  xi_2] 
[  phi_3  xi_3] 
[  phi_4  xi_4] 


phi.i  is  the  delay,  and  xi_i  is  the  channel  gain  for  path  i 
mxi:  the  largest  delay 

W  :  2\beta\omega\sinc(2\pi\omega\delta  T  n)  , 

kk  :  a  constant  that  allow  us  to  manipulate  the  eigenvalues  of  th 
kernel 


V  =  zeros (size (r .alpha) ) ; 

Void  =  V; 

b  =  randn (size (r .alpha)  )  ; 
b  =  b/norm(b) ; 
r  =  b; 
p  =  b; 
flag  =  l; 
th  =  le-5 ; 
arg  =  1; 
i  =  0; 

while  (i<30)  &  (arg  >  le-5) 
i  =  i  +  l; 
k=0; 

while  flag, 
k=k+l; 

z  =  Rxb(r. alpha, s,xi,mxi,W,kk,p) ; 

v  =  (r ,*r)/(p,*z) ; 

V  =  V  +  v*p; 
rnew  =  r  -  v*z; 

f printf1 T> \b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b > ) ; 
f pr intf ( ’ \b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\b\bMineig_cg_method :  ) ; 

f printf  (’7.d  7.d  7.f  ’  ,i,k,temp) ; 
if  temp  <  th 
flag  =  0; 

end 

mu  =  (rnew ’ *rnew) / (r ’ *r ) ; 
r  =  rnew; 
p  =  r  +  mu*p; 
end 

V  =  V/norm(V) ; 
arg  =  norm(V-Vold) ; 

Void  =  V; 
th  =  th/1.25; 
r  =  V; 
p  =  V; 
flag  -  i; 
end 

y  =  Rxb (r. alpha, s,xi, mxi, W, kk, V) ; 

ev  =  norm(y)-kk; 
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Land-Mobile  Model  Autocorrelation  Function 

function  r  =  Cal_r (beta,Bd,n,deltaT) 

•/,  r  =  Cal_r(beta,Bd,n, delta!) 

V.  fading  channel  autocorrelation 

•/,  This  program  calculates  the  land-mobile  fading  cn 

•/,  coeffs. 

•/. 

I  Run  the  composite  mapping  algorithm  to  make  a  positive  sequence 
•/,  Modified  from  Dr  Moon’s  COMPMAP.m 


n=n/2; 

N  =  4*n; 


q  =  n-1;  t  =  (0:q)’; 


arg 
h  = 

X  = 


:  2*pi*Bd*deltaT 
jesselj (0,arg*t) ; 

[h:  zeros (N-2*q-l > 1) 


h(q+l:-l:2)3 


•/,  Conjugate  even  extension 


'/,  Run  the  composite-mapping  algorithm 
converged  *  0;  nnmit.r  -  0;  makiter  ■  2000, 
while  (numiter  <  maxiter) 
numiter  =  numiter+1 , 

X  =  real(fft(x))  ; 

idx  =  X  <  °; 

X(idx)  -  zeros (size(X(idx) )) ; 
newx  =  real (if f t (X) ) ; 
newx(q+2:N-q)  =  zeros (M-2*q-l,l) ; 
if (norm(x-newx)  <  l.e-10) 


•/,  Enforce  Property  2:  positive  sequence 

*/,  Enforce  Property  1:  length  2q+l 
y,  check  for  convergence 


break; 

end 

x  =  newx; 
end 

h=x(l :  (q+D)  *> 
h=h/h(l) ; 
a  =  zeros (n,l) ; 
r=  [h ;  a]  ; 
r=r* (1-beta) ; 


Frequency  Containment 

function  h  =  Cal_W (bet a , omega , n , delt aT) 

*/,  h  =  Cal _W (beta, omega .n.deltaT) 

l  Thin  program  calculates  the  conffs  of  sine  function  for  freq 
'/,  containment . 

\  h  ,  2\beta\omega\sinc(2\pi\omega\delta  T  n) 

\  Rm  the  composite  mapping  algorithm  to  make  a  positive  sequence 
'/,  Modified  from  Dr  Moon’s  COMPMAP.m 
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N  =  4*n;  q  =  n-1;  t  =  (0:q)’, 

argl  =  2*beta*omega; 

arg2  =  6.28319*omega*deltaT; 

h  =  sinc(arg2*t) ;  , 

x  =  [h;  zeros (N-2*q-l, 1) 5  h(q+l:-l:2; J ; 


%  Conjugate  even  extension 


•/,  Run  the  composite-mapping  algorithm 
converged  =  0;  numiter  =  0;  maxiter 
while  (numiter  <  maxiter) 
numiter  =  numiter+1 ; 

X  =  real(fft(x)) ; 


2000; 


idx  =  X  <  0; 

X(idx)  =  zeros (size (X(idx) )) i 

newx  =  real (if f t (X) ) ; 
newx(q+2:N-q)  =  zeros (N-2*q-l, D > 
if (norm (x-newx)  <  1 . e- 10) 


'/,  Enforce  Property  2:  positive  sequence 

*/.  Enforce  Property  1:  length  2q+l 
V,  check  for  convergence 


break; 


end 

x  =  newx; 
end 

h=x( 1: (q+1)) J 
h  =  h/h(l) ; 
h  =  argl*h; 


Synchronization  Delay  Distribution 


function  S  -  SPmu(s,K,mu) 

'/,  function  S  =  SPmu(s,mu) 


•/. 

•/. 

•/. 


i.  program  is  used  to  calculate  the  eMective  sign.l/f liter 
-or dine  to  synchronization  delay  distribution 


'/,  mu:  0->delta,  l->uniform,  2->Gaussian 
•/,  K:  span  of  the  random  variable 

y# 

*/,  S  [n]  =  sum^k  s(n-k)pmu(k) 


if  mu  ==  0 
S  =  s; 
else 

N  =  size (s , 1) ; 

S  =  zeros (N,l) ; 

M  =  (K-l)/2; 

A  =  zeros (M,l) ; 

Sa  =  [A ; s ; A] ; 
if  mu  ==  1 

pmu  =  ones (K, 1) ; 
else 

k  pmu  =’ gaussian.pdf ( (7 . 38/K) * (k-M-1) ) ; 
end 
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pmu  =  pmu/sum(pmu) ; 
for  n=l:N 
for  k  =  1:K 

S(n)  =  S(n)+Sa(n+K-k)*pmu(k); 
end 
end 

end 


Gaussian  Density  Function 

function  pdf  =  gauss ian_pdf (x) 

V,  pdf  =  gaussian_pdf  (x) 

V,  This  function  returns  the  value  of  Gaussian  distribution 
*/,  with  zero  mean  and  unit  variance 
•/. 


[r ,  c]  =  size  (x) ; 

s  =  r  *  c; 

x  =  reshape(x,  i,  s) ; 
pdf  =  zeros (1,  s) ; 
k  =  f ind(isnan(x)) ; 
if  (any(k)) 

pdf(k)  =  NaN  *  ones(l,  length(k)); 
end 

k  =  find  ("isinf(x)); 

if  (any(k))  ,>  .,/,v 

pdf (k)  =  (2  *  pi)~("  1/2)  *  exp(-  x(k)  .  2  /  )> 

end 

pdf  =  reshape (pdf,  r,  c) ; 


Vector  Normalization 

function  [V]=normal_v(W,deltaT) 

V,  [V]=normal_v(W) 

v 

y  Both  V,  and  W  are  (column  or  row)  sample  vectors  of  a 
%  continuous  function,  where  V  is  an  normalized  version 
y,  of  W  such  that  |V I  =  1  (in  continuous  sense) 

•/. 


D  =  size(W); 
if  D(l)  ==  1 

V  =  W. /sqrt (W*W5 *deltaT) ; 

elseif  D(2)  ==1 

V  =  W . /sqrt (W* *W*deltaT) ; 

else 

V  =  W; 
end 
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C  Codes 


Matrix- Vector  Multiplication 


/ 


************************************** 
V  =  Rxb(r_alpha,s,xi,mxi,W,k,b) 


************************************* 


Do  the  matrix_vector  multiplication: 
[(l_beta)R_ alpha  o  Psi(\tau,\lambda) 
=  R*b 


+  beta*w*sinc [w(\tau-\lambda)]+kl] *b 


We  have: 

R_alpha(m,n)  =  r_alpha( |m-n| ) 
xi  =  [  phi_l  xi_l] 

[  phi_2  xi_2] 

[  phi_3  xi_3] 

[  phi_4  xi_4] 

phi.i  is  the  delay,  ahd  ti.i  is  the  charnel  gaia  lor  path  i 
mxi  is  the  largest  delay 

Psi(m,n)  =  \sum_i  lxi.il‘2  s (m-\phi_i) s (n-\phi) 

V(n)  =  w*sinc[w(n)] 


note: 


>  m 


n 


V(n)  =  \sum_m  R(m,n)*b(m) 

=  \sum  m[(i_beta)r_alpha(|m-n| )* 

{\sum_i  I  xi.il ‘2  s (m-\phi_i) s (n-\phi)  } 
+  bet a*W ( I m-n I ) +k*delta (m , n) ] *b (m) 


bet.  Bd  1-beta  are  absorbed  ih  r.alpha  and  V  respectively  before 
this  function  is  called 

. . — ............... — .......**«***♦* — **'*' 

# include  <math.h> 

♦include  "mex.h" 


/*  Input  Arguments  */ 

#def ine  ir_alph  prhs [0] 
♦define  is  prhs[l] 


#def ine  ix  prhs[2] 
#def ine  imxi  prhs[3] 
#def ine  iW  prhs[4] 
#def ine  ik  prhs[5] 
#def ine  ib  prhs[6] 


/*  Output  Arguments  */ 


#define  ov  plhs[0] 


void  mexFunctionC  int  nlhs,  mxArray  *Plhs[], 
int  nrhs,  const  mxArray*prhs  □  ) 


double  *r_alpha, *s , *xi , *xil , *mxi 

int  i,n,m,M,N,I; 

int  sshift,  mphi,  shift; 


*W ,  *k ,  *b ,  *v , *vn , *S , *SS ; 


double  R,  temp; 
mxArray  *Sarray; 


r .alpha  =  mxGetPr (ir .alpha) ; 
s  =  mxGetPr (is) ; 
xi  =  mxGetPr (ixi) ; 
mxi  =  mxGetPr (imxi) ; 

W=  mxGetPr (iW) ; 
k  =  mxGetPr (ik) ; 
b  =  mxGetPr (ib) ; 

M  =  mxGetM(ib) ; 

ov  =  mxCreateDoubleMatrix (N , 1 ,mxREAL) ; 
v  =  mxGetPr (ov); 
f  or (n=0 ; n<N ; n++)  *(v+n)=0; 

I  =  mxGetM(ixi) ; 
xil  =  xi+I; 

sshift  =  0; 

/*  Prepare  for  allowing  shifts  in  s  for  path  delays  */ 

if(I!=lH 

ii(*mxi  >  OH 
sshift  =  *mxi; 

shift  =  N+sshift;  __,rN 

Sarray  =  mxCreateDoubleMatrix(shif t , 1 ,mxREAL) ; 

S  =  mxGetPr (Sarray) ; 

for (n=0 ; n<sshif t ; n++)  *(S+n)  =0; 

for (m=0 ; n<shif t ;n++ ,m++)  *(S+n)  =  *(s+m); 

> 

else{ 

shift  =  N  -  *mxi; 

Sarray  =  mxCreateDoubleMatrix(shift,i,mxREAL) , 

S  =  mxGetPr (Sarray) ; 


for(n=0;n<N;n++)  *(S+n)  -  *(s+n); 
for(;n<shift;n++)  *(S+n)  =  0; 

> 

> 

else{ 

Sarray  —  mxCreateDoubleMatrix(N,l,mxREAL) , 

S  =  mxGetPr (Sarray) ; 

for (n=0 ; n<N ; n++)  *(S+n)  =  *(s+n); 

> 

/*  Divided  into  three  parts  to  avoid  abs(m-n)  */ 
for  (n=0 ;  n<N ;  n++)  { 
vn  =  v+n; 
for  (m=0;m<n;m++)-( 

R  =  0; 

for(i=0;i<I;i++H 

shift  =  s shift  -  *(xi+i); 

SS  =  S+shift; 

temp  =  * (xil+i) ** (SS+m) ; 

temp  =  temp**(SS+n) ; 

R  =  R  +  temp; 

> 

shift  =  n  -  m; 

R  =  R**(r_alpha+shift) ; 

R  =  R  +  *(W+shift); 

R  =  R**(b+m) ; 

*vn  =  *vn  +  R; 

> 


for  (n=0;n<N;n++M 
vn  =  v  +  n; 

R  =  0; 

for(i=0;i<I;i++H 

shift  =  n  +  sshift  -  *(xi+i); 
SS  =  S+shift; 
temp  =  *SS**SS j } 
temp  =  temp** (xil+i) ; 

R  =  R  +  temp; 

> 

R  =  R** (r_alpha) ; 

R  =  R  +  *(W) ; 

R  =  R  +*k; 
r  =  R**(b+n) ; 

*vn  =  *vn  +  R; 

> 

for (n=0 ; n<N ; n++) { 
vn  =  v+n; 

for  (m=n+ 1 ;  m<M ;  m++)  ■( 

R  =  0; 
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for(i=0;i<I;i++){ 

shift  =  sshift  -  *(xi+i); 
SS  =  S+shift; 
temp  =  *(xil+i) **(SS+m) ; 
temp  =  temp**(SS+n) ; 

R  =  R  +  temp; 

} 

shift  =  m  -  n; 

R  =  R**(r_alpha+shift) ; 

R  =  R  +  *(W+shift) ; 

R  =  R**(b+m) ; 

*vn  =  *vn  +  R; 

} 

} 

mxDestroyArray (Sarray) ; 
return; 

> 


